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ABSTRACT 

We performed a detailed analysis of elemental abundances, dust features, and polycyclic aro¬ 
matic hydrocarbons (PAHs) in the C-rich planetary nebula (PN) Wray 16-423 in the Sagittarius 
dwarf spheroidal galaxy, based on a unique dataset taken from the Subaru/HDS, MPG/ESO 
FEROS, //5T/WFPC2, and Spitzer/ IRS. We performed the first measurements of Kr, Fe, 
and recombination O abundance in this PN. The extremely small [Fe/H] implies that most 
Fe atoms are in the solid phase, considering into account the abundance of [Ar/H]. The 
Spitzer /IRS spectrum displays broad 16-24 pm and 30 pm features, as well as PAH bands 
at 6-9 pm and 10-14/im. The unidentified broad 16-24/im feature may not be related to iron 
sulfide (FeS), amorphous silicate, or PAHs. Using the spectral energy distribution model, 
we derived the luminosity and effective temperature of the central star, and the gas and dust 
masses. The observed elemental abundances and derived gas mass are in good agreement with 
asymptotic giant branch nucleosynthesis models for an initial mass of 1.90M o and a metal- 
licity of Z-0.004. We infer that respectively about 80%, 50%, and 90% of the Mg, S, and 
Fe atoms are in the solid phase. We also assessed the maximum possible magnesium sulfide 
(MgS) and iron-rich sulfide (Fe50S) masses and tested whether these species can produce the 
band flux of the observed 30 pm feature. Depending on what fraction of the sulfur is in sulfide 
molecules such as CS, we conclude that MgS and Fe50S could be possible carriers of the 
30 pm feature in this PN. 

Key words: ISM: planetary nebulae: individual (Wray 16-423) — ISM: abundances — ISM: 
dust, extinction 


1 INTRODUCTION 

Planetary nebulae (PNe) represent the final stage in the evolution 
of initially ~1-8M 0 stars. During their evolution, these mass stars 
eject a large amount of their mass into circumstellar shells via 
strong stellar winds. Eventually, the ejecta blend with interstel¬ 
lar media (ISM). The evolution history of the progenitors is im¬ 
printed in the central stars of PNe (CSPNe), in the atomic gas, 
dust, and molecule-rich circumstellar shells. Using the emission 
lines of gas and thermal radiation from dust grains, one can si¬ 
multaneously obtain abundances of elements and dust, as well as 
the ejected gas/dust masses from the PN progenitors. An investi¬ 
gation of these materials provides insight into stellar evolution and 
the ISM material-recycling mechanism of host galaxies. Thus, PNe 


* Based on data collected by the Subaru Telescope, operated by the Na¬ 
tional Astronomical Observatory of Japan (NAOJ) under program IDs: 
S12A-126S and S14A-174 (PI: M. Otsuka). 

f Based on observations made with the MPG/ESO 2.2-m Telescope at the 
La Silla Observatory under program ID 91.D-0055A (PI: M. Otsuka). 
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provide a unique laboratory for the study of stellar and galaxy evo¬ 
lution. 

Elemental abundances, the luminosities of CSPNe, and gas 
and dust masses are essential to understanding stellar and galaxy 
evolution; the latter three parameters are more accurately ob¬ 
tained through observations of extragalactic PNe, rather than 
Galactic PNe, because the distance is well determined. Currently, 
there are over 4,200 known extragalactic PNe (from Table 1 of 
Reid 2012). In the galaxy that is nearest to us (24.8±0.8 kpc; 
Kunder & Chaboyer 2009), the Sagittarius dwarf spheroidal galaxy 
(Sgr dSph), four PNe have been confirmed to date; BoBnl, StWr2- 
21, Hen2-436, and Wray 16-423 (Walsh et al. 1997; Dudziak et al. 
2000; Zijlstra et al. 2006a; Kniazev et al. 2008; Otsuka et al. 2008, 
2010 , 2011 ). 

Sgr dSph PNe offer obvious advantages over other extragalac¬ 
tic PNe in terms of the study of stellar evolution, such as those 
found in Magellanic Clouds (MCs); Sgr dSph PNe are closer to us 
and therefore appear brighter than other extragalactic PNe, allow¬ 
ing the detection of rare elements such as neutron (ra)-capture ele¬ 
ments. The final yields of neon and n-capture elements are largely 
concerned with the 13 C pocket region that forms in the He-rich in- 
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tershell during the thermal-pulse asymptotic giant branch (AGB) 
phase. The extent of the 13 C pocket is an important parameter in 
AGB nucleosynthesis models; however, it remains unconstrained. 
An assessment of n-capture and neon abundances, and the progen¬ 
itor mass derived from Sgr dSph PNe could be used to test the 13 C 
pocket mass adopted in AGB models. 

The C/O ratio is an important parameter with regard to CSPN 
and nebula characterisation and initial mass constraints. The C 
and O abundances using recombination lines (RLs) are sometimes 
larger than those using collisionally excited lines (CELs; e.g., see 
Liu 2006, for details). Moreover, the discrepancy factor between 
the RL and CEL abundances is generally different for O and for C 
(e.g., see Wang & Liu 2007). As a consequence, the C/O ratio de¬ 
termined from dividing the RL C abundance by the CEL O abun¬ 
dance might not be indicative of the actual C/O ratio. Therefore, 
we should at the very least determine the C/O ratios using the same 
type of emission lines. According to Delgado-Inglada & Rodriguez 
(2014), the uncertainty of the RL C/O ratio (±0.06 dex) is con¬ 
siderably smaller than that of the CEL C/O ratio (±0.2 dex). An 
important factor of the CEL C/O uncertainty is the normalization 
between UV and optical fluxes since CEL C abundances are deter¬ 
mined using the UV [Cm,IT 1906/09 A lines whereas we measure 
the CEL O abundances using the optical [O in] TT4959/5007 A and 
[O n] TT3726/29 A. By using C and O RLs in the optical, we do not 
have to worry about this difficult flux normalization and one can 
simultaneously determine C and O abundances and then the C/O 
ratio, as performed in Sgr dSph PNe BoBnl (Otsuka et al. 2010) 
and Hen2-436 (Otsuka et al. 2011). 

The Hubble Space Telescope ( HST ) can resolve the central 
stars of Sgr dSph PNe, as reported in Zijlstra et al. (2006a). There¬ 
fore, the luminosities of the CSPNe can be measured directly, with¬ 
out the need for empirical methods. The mass of the progenitor star 
can be estimated by plotting the measured luminosity and temper¬ 
ature of the CSPN on theoretical post-AGB evolution tracks. Thus, 
using this approach, the observed elemental abundances and gas 
mass can be compared to predicted values from AGB nucleosyn¬ 
thesis models. 

Previously, we carried out multiwavelength observational 
studies of the Sgr dSph PNe BoBnl (Otsuka et al. 2010) and 
Hen2-436 (Otsuka et al. 2011). While all Sgr dSph PNe are 
metal-poor (Walsh et al. 1997; Dudziaketal. 2000; Zijlstra et al. 
2006a), BoBnl is extremely metal-poor when comparing to Hen2- 
436 or Wrayl6-243 (Otsuka et al. 2015; Zijlstra et al. 2006a; 
Kniazev et al. 2008). Our study on chemical abundances and dust in 
these Sgr dSph PNe has provided us with much information on PNe 
and their progenitor’s evolution. Otsuka et al. (2010) measured the 
rare n-capture elements Xe and Ba abundances for first time in ex- 
tragalactic PNe. Based on the abundance pattern of BoBnl, they 
proposed that the progenitor might be a binary. Otsuka et al. (2011) 
measured the n-capture Kr abundance for the first time in Hen2- 
436 and found that the abundance pattern is very similar to that 
of Wray 16-423. Both studies furthermore also derived dust masses 
and discussed the dust production in metal-poor environments. In 
this paper, we focus on the Sgr dSph PN Wray 16-423. 

Wray 16-423 was discovered near the core of the Sgr dSph by 
Zijlstra & Walsh (1996), based on direction and radial velocities. 
Ibataetal. (1995) reported an average heliocentric radial veloc¬ 
ity for the Sgr dSph of 140±2km s _1 , with an intrinsic velocity 
dispersion of 11.4±0.7km s ^ 1 (Ibataetal. 1997). The radial ve¬ 
locity of Wrayl6-423 was determined to be +133.1±2km s~' by 
Zijlstra & Walsh (1996). The elemental abundances of the nebula 
and the luminosity and effective temperature of the CSPN were 


calculated by Walsh et al. (1997); these values were refined by the 
photoionisation models of Dudziak et al. (2000), Gesicki & Zijlstra 
(2003), and Zijlstra et al. (2006a, see Table 4 and Section 3.4 for 
a comparison of these values). Dudziak et al. (2000) calculated 
the respective C and O abundances from the C RLs and the O 
CELs. The C/O ratio using the same type of emission lines is 
still unknown. Dudziak et al. (2000) expected an [Fe/H] abundance 
of —0.5, although they did not report the detection of Fe lines. 
Zijlstra et al. (2006a) determined the magnitude of the CSPN to be 
niv=18.85±0.20, based on /LST/F547M measurements. However, 
this estimate seems to be too bright. Indeed, assuming a spectral 
energy distribution defined by an m v = 18.85 CSPN and adopting 
an effective temperature of 107 000 K (Zijlstra et al. 2006a), the lu¬ 
minosity L, should be ~12 600L o . However, this does not match 
their estimated value of L»=4350L G obtained using a photoionisa¬ 
tion model. This discrepancy and a revision of the CSPN magnitude 
is discussed in Section 2.4. 


Stanghellini et al. (2012) reported that the Spitzer mid- 
infrared spectrum of Wrayl6-423 displays carbonous dust features. 
As we report later, the broad 30 pm feature is seen in this PN. 
The broad 30 Tm feature has been found in many C-rich objects. 
However, the carrier of the broad 30 pm feature is unidentified; 
Hony et al. (2002a) proposed that the feature originates from mag¬ 
nesium sulfide (MgS) whose shapes can be described by a con¬ 
tinuous distribution of ellipsoids (CDE). Since then, CDE MgS 
grains have been the major candidate for the 30 pm feature (e.g., 
Zijlstra et al. 2006b; Sloan et al. 2014). However, the recent study 
of Zhang et al. (2009) casts doubt on this identification for the 
30pm feature. They demonstrated that the MgS mass formed from 
the available S atoms, could not account for the observed strength 
of the 30 pm feature in the Galactic proto-PN HD56126, even if 
all the dust grains existed as MgS. CDE iron-rich sulphide, such 
as Mgo. 5 Feo. 5 S (Fe50S), also show a broad feature around 30/rm, 
similar to the CDE MgS feature (Messenger et al. 2013). Using the 
data of Wray 16-423, we can test whether the strength of the 30yum 
feature can be explained by thermal emission from the available 
MgS or Fe50S grains. Since the distance to Wray 16-423 is well de¬ 
termined, one can directly determine the luminosity of the CSPN as 
the gas and dust heating source and set the size of the dusty nebula. 
Therefore, the solid phase Mg, S, and Fe masses and the MgS and 
Fe50S grain masses would be more credible by using the accurately 
measured Mg, S and Fe abundances. 


For the above reasons, we performed a detailed spectroscopic 
analysis of Wray 16-423 based on multiwavelength data taken from 
Subaru/HDS, MPG ESO/FEROS, Spitzer/ IRS, and W57/WFPC2: 
( 1 ) to determine and refine nebular elemental abundances, in par¬ 
ticular, the C and O recombination abundances and Fe and slow 
n-capture process (j-process) elements; ( 2 ) to investigate dust fea¬ 
tures and PAH bands and compare these among LMC PNe; (3) to 
derive the physical conditions of the CSPN, gas, and dust by con¬ 
structing the spectral energy distribution (SED) model using the 
radiative transfer code Cloudy (Ferland et al. 1998); (4) to estimate 
the progenitor mass by comparing the observed elemental abun¬ 
dances and estimated gas mass with AGB model predictions; and 
(5) test how much masses of the MgS and Fe50S grains are neces¬ 
sary to reproduce the broad 30 pm feature. 
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Table 1. Summary of optical spectroscopic observations. 

Date Telescope/Instrument Wavelength Exposure time 

2012/07/04 Subaru 8.2-m/HDS 4300-7100A 2x1800s, 360s 

2013/06/18 MPG/ESO 2.2-m/FEROS 3500-9200 A 3x3000 s, 2x300 s 
2014/07/09 Subaru 8.2-m/HDS 3680-5400 A 2x1800 s, 180 s, 30 s 

2 OBSERVATIONS AND DATA REDUCTIONS 

2.1 Subaru 8.2-m Telescope HDS spectroscopy 

Optical high-dispersion spectra were taken using the High- 
Dispersion Spectrograph (HDS; Noguchi et al. 2002) attached to 
the Nasmyth focus of the 8.2-m Subaru Telescope on 2012 July 4 
(Prop.ID: S12A-126S, PI: M. Otsuka) and 2014 July 9 (Prop.ID: 
S14A-174, the same PI). HDS observations are summarised in Ta¬ 
ble 1. 

The weather conditions during the exposure for both runs were 
stable, and the visibility was ~0.5" (2012) and ~0.8" (2014), mea¬ 
sured using a guider CCD. An atmospheric dispersion corrector 
(ADC) was used to minimise the differential atmospheric disper¬ 
sion through the broad wavelength region. The slit width and length 
were set to 1.2" and 6", respectively, and the 2x2 on-chip binning 
mode was selected. The resolving power (A/AA) reached 33 500, 
which was measured from the average full width at half maximum 
(FWHM) of over 600 narrow Th-Ar comparison lines obtained for 
wavelength calibrations. For the flux calibration, blaze function 
correction, and airmass correction, we observed the standard star 
HD 184597. 

We reduced the data with the echelle spectra reduction pack¬ 
age ECHELLE and the two-dimensional spectra reduction package 
TWODSPEC in IRAF 1 . To keep a good signal-to-noise (S/N) over 
the observed wavelength range, we extracted a region of 15 pixels 
along the slit (corresponding to 4.14") and summed all fluxes. For 
the flux calibration, we reduced the ESO archival spectra (~0.3- 
2.5 pm) of HD 184597, taken using XSHOOTER (Vernet et al. 
2011, Prop.ID: 60.A-9022C). We corrected the reduced 0.3-2.5 gm 
spectrum of HD 184597 to match the B- and V-band magnitudes re¬ 
ported in H0g et al. (2000, 6=6.81 and V=6.88), and then we con¬ 
structed a wavelength versus A6-magnitude table in small wave¬ 
length increments (0.4 A, for this study) for performing flux cali¬ 
bration of the Wray 16-423 spectra. 

The resulting S/N ratio in the sky-subtracted object frame of 
Wray 16-423 is >5 for the continuum and >8 at the intensity peaks 
of the emission lines. 

The resultant HDS spectrum is presented in Figures Al and 
A2 of Appendix A. 

2.2 MPG ESO 2.2-m Telescope FEROS spectroscopy 

Optical high-dispersion spectra (3500-9200 A) were obtained us¬ 
ing the Fiber-fed Extended Range Optical Spectrograph (FEROS; 
Kaufer et al. 1999) attached to the MPG/ESO 2.2-m Telescope, La 
Silla, Chile on 2013 June 18 (Prop.ID: 91.D-0055A, PI: M. Ot¬ 
suka), as a follow-up up to HDS observations; auroral [On] lines, 
nebular [S m], [Cl iv], and [Ar m] lines, and the Paschen discontinu¬ 
ity were measured, and line identifications and flux measurements 

1 IRAF is distributed by the National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy 
(AURA), Inc., under a cooperative agreement with the National Science 
Foundation. 



Wavelength (jum) 


Figure 1. Spitzer /IRS spectrum of Wrayl6-423. As-observed atomic lines, 
PAHs, and broadband features are indicated. 

were cross-checked for both HDS and FEROS spectra. The FEROS 
observation is summarised in Table 1. 

The weather conditions were stable and clear throughout the 
night, and the visibility was ~ 1.1" measured from the guider CCD. 
FEROS’s fibers use 2.0" apertures, to provide simultaneous spectra 
of the object and sky frames. The ADC was used during the ob¬ 
servation. We selected a 2x2 on-chip binning mode and low gain 
mode 2 . The resolving power reached ~37 200 measured from the 
average FWHM of over 100 Th-Ar comparison lines. We also took 
5x3600 sec dark frames to correct for the zero intensity of longer 
exposure-time frames. For the same reason as described in the HDS 
spectroscopy, we summed over the ±3 pixels centered on the fibre 
position. We confirmed that the extracted spectra detected ~65 % 
of the light falling into the 2" aperture by comparing to the spec¬ 
tra made by summing up the full aperture. For flux calibration, we 
observed the standard star HR7950, before and after the Wray 16- 
423 exposure. We referred to the HR7950 spectrum of Valdes et al. 
(2004) for constructing a wavelength versus A6-magnitude (0.4 A 
step). 

We reduced the data using the same method as that adopted for 
the HDS data. The resulting S/N ratio in the sky-subtracted object 
frame is >2 on the continuum and >5 at the intensity peaks of the 
emission lines. 

2.3 Spitzer/IRS spectroscopy 

We analysed the archival mid-infrared spectra taken by the 
Spitzer/Infrared Spectrograph (IRS; Houck et al. 2004) with the 
SL (5.2-14.5 pm, the slit dimension: 3.6"x57") and the LL mod¬ 
ules (13.9-39.9yum, 10.6"xl68"). The data were originally taken 
by L. Stanghellini (AOR Key: 25833216) on 2008 June 9. We 
processed them using the reduction package SMART v.8.2.9 
(Higdon et al. 2004) and IRSCLEAN v.2.1.1 provided by the 
Spitzer Science Center. We performed additional correction of the 
IRS spectrum by scaling it up to match the flux density at the band 
W4 (d c =22.09 pm, AA=A.Wpm) of the Wide-field Infrared Sur¬ 
vey Explorer ( WISE', Wright et al. 2010), using a constant factor 
of 1.069. 

The resulting spectrum is presented in Figure 1. The spectrum 

2 The gain and readout-noise in the case of 2x2 binning and low gain 
mode were measured to be 4.03 e ADU and 5.68 e using the IRAF task 
FINDGAIN. 
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Off-Set (HA) [arcsec] 

Figure 2. HST/WPVC2 F656N (Ha) image of Wray 16-423. The levels of 
the contour plots are from 3-cr above the mean background intensity to the 
peak intensity of the object. The radius of the green dashed circle is 1.3". 


Table 2. Photometry of the HST /WFPC2 images. 


Band 

A c 

(A) 

Portion 

Fa 

(erg s -1 cm -2 A -1 ) 

F547M(V) 

5483.86 

CSPN only 

(3.75±0.34)xl0 -17 



whole nebula 

(1.1 l±0.04)x 10 15 

F656N(Ha) 

6563.76 

whole nebula 

(9.61±0.15)xl0 -14 


Note - The V-band magnitude of the CSPN is 19.98±0.20, determined 
from the F547M image. 

clearly shows the 6-9gm and 11.3/tm PAHs, the broad 11/jm, lb- 
24 pm, and 30 /tm bands. 


2.4 HST/WFPC2 photometry and the Ha/H/J fluxes 

We reduced the HST /Wide Field Planetary Camera 2 (WFPC2) 
F656N (Her, T C =6563.76A, Ad=53.77 A) and F547M images 
(V-band, /l c =5483.86 A, AT=205.52A) using the MultiDrizzle 
(Koekemoer et al. 2003) on PYRAF. The data were taken on 2002 
July 25 (Prop.ID: 9356, PI: A. Zijlstra). We set the plate scales to a 
constant 0.025" pixel -1 . The resultant F656N image is displayed in 
Figure 2. The Wray 16-423 image shows an inner cylindrical struc¬ 
ture, surrounded by an elliptical nebula shell that extends radially 
to —1.3", as indicated by the green dashed circle. In the F656N and 
F547M images, we measured the count rates (cts) within the aper¬ 
ture radius of 1.3", corresponding to the sum of the CSPN's and 
nebula’s cts. We measured average (background) cts in an annulus 
with inner radius 2.1" and outer radius 2.4", centered on the CSPN. 
Finally, we converted the cts into flux density, using the WFPC2 
photometric zero-points in erg s -1 cm -2 A -1 cts -1 . The results are 
summarised in Table 2. 

To measure the pure Ha flux using the F656N flux density, one 
must remove the contribution from both the local continuum and 
the [N ii] A 6548 A. We used the FEROS spectrum to estimate these 
contributions. Taking into account the F656N filter transmission 



Radial Velocity (km s -1 ) 


Figure 3. (from left to right) The emission line profiles of the 
[O n] A 3726.03 A, H y 4340.46 A, and [O hi] A 4958.91 A, from HDS obser¬ 
vations. The wavelength is converted into radial velocity with respect to the 
systemic radial velocity of +133.12 km s -1 . The flux density is normalised 
to the flux density at the peak. 


characteristics, we compared the F A (HST, F656N) with the coun¬ 
terpart FEROS spectrum, i.e., F A (FEROS,F656N). Finally, we de¬ 
termined the scaling factor F A (HST,F656N)/F A (FEROS,F656N) = 
2.383 3 . 

Using this scaled FEROS spectrum, we measured a pure 
F(Ha) of (4.70±0.03)xl0 -12 erg s -1 cm -2 and a pure F(H0) of 
(1.53±0.02)xl0 -12 erg s -1 cm -2 . We used the measured pure _F(H/3) 
to normalise the flux in the flux density scaled Spitzer /IRS spec¬ 
trum, as explained in the previous section. 

Next, we measured the V-band flux density of the CSPN in 
the F547M image. In the flux density measurements of the CSPN, 
we need to select the background area more carefully because the 
CSPN is embedded in a bright nebula and therefore the CSPN 
magnitude largely depends on the background region adopted. 
The background cts adopted in the flux density measurement of 
the whole nebula are not suitable for the flux density measure¬ 
ments of the CSPN. For this reason, we measured cts in four po¬ 
sitions near (~0.2-0.3" away) the CSPN and used the average cts 
as the sky background. The measured flux density of the CSPN 
in the F547M band corresponds to the V-band magnitude my of 
19.98±0.20, which is about 1 magnitude fainter than Zijlstra et al. 
(2006a, 77iy=18.85±0.20). The difference would be attributed to 
the adopted background region; the my is 18.73±0.10 when as the 
background we use the average cts measured in the area of the an¬ 
nulus centered on the PN with inner and outer radii of 2.1" and 
2.4", respectively. 


3 RESULTS 

3.1 Reddening correction, line-flux measurements, and 
radial velocity 

The measured line fluxes in the obtained spectra were de-reddened 
using the following formula: 

1(A) = F(A) X IQ c0mi+M) 9 (i) 


3 We would estimate the Wray 16-423’s nebula to be >2.8" at the 1.1" 
seeing, assuming that its actual nebula size is 2.6". From our report in Sec¬ 
tion 2.2, our FEROS spectra detected ~0.65x(2.0"/2.8")x 100~46% or 
less of the light from Wray 16-423. In the actual observation, the seeing was 
not a constant 1.1" and also the thin clouds were passing during the expo¬ 
sures. This is the reason why we need such a correction factor. 
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where 1(A) is the de-reddened line-flux, F(A) is the observed line- 
flux, f(A) is the interstellar extinction function at A computed by 
the reddening law of Cardelli et al. (1989) with R v = 3.1, and c(H6) 
is the reddening coefficient at H ft. For the Spitzer/IRS spectrum, we 
adopted the reddening law of Fluks et al. (1994). 

We measured c(H/3) by comparing the observed Balmer 
line ratios of Fly and Ho- to H/3 with the theoretical ratios of 
Storey & Hummer (1995) for an electron temperature r E =10 4 K 
and electron density /2 f =10 4 citT 3 , under the Case B assumption. 
The Her line in the FIDS short-exposure spectrum taken in the 
2012 run was affected by cosmic rays. Therefore, we used a c(H/3) 
value of 0.132±0.035, determined from the F(Hy)/F(H/3) ratio. 
This c(H/3) value is in excellent agreement with that in the 2014 
HDS run (0.133±0.063). In the FEROS spectrum, we obtained a 
c(H/3) value of 0.110±0.022 from the F(H»)/F(H/?) ratio. 

We fitted the emission lines in the HDS and FEROS spectra 
using multiple Gaussian components. Figure 3 shows the emis¬ 
sion line profiles of the [On] A 3726.03 A, Hy A 4340.46 A, and 
[O hi] A 4958.91 A. We obtained a heliocentric radial velocity of 
+ 133.12±0.28 km s _1 (RMS of the residuals: 3.17 km s~‘) from the 
central wavelengths of all the lines detected in the HDS and FEROS 
spectra. Most line profiles can be fitted by a single Gaussian com¬ 
ponent, whereas several lines show double peaks or blue-shifted 
asymmetry, as seen in [On] and [Om], 

We list the observed wavelengths and de-reddened relative 
fluxes of each Gaussian component, indicated by the identification 
(ID) number in Table Al of Appendix A, with respect to the de- 
reddened H/l flux /(Hj8) of 100. For lines composed of multiple 
components, we list the de-reddened relative fluxes of each com¬ 
ponent, as well as the sum of these components (indicated by ’T’). 
In the ions showing multiple intensity peaks such as [O ii] 43726 A, 
one could derive electron temperatures and densities in each Gaus¬ 
sian component. In our data, it would be difficult to accurately de¬ 
rive their ionic abundances relative to the H + in each component 
because each velocity component of the detected lines are differ¬ 
ent. Since our aim is to investigate average ionic and elemental 
abundances in the nebula, we used the integrated fluxes of each 
Gaussian component. We list the results of the flux measurements 
for the Spitzer/IRS spectrum in Table A2 of Appendix A, where we 
adopted the c(H/3) of 0.110±0.022, which is the same value applied 
for the FEROS spectrum. 

3.2 Plasma diagnostics 

In the following line-diagnostics and subsequent ionic element cal¬ 
culations, the adopted transition probabilities, effective collision 
strengths, and recombination coefficients are the same as those 
listed in Tables 7 and 11 of Otsuka et al. (2010). 

We calculated the electron densities n £ and temperatures T £ 
using the diagnostic line ratios of RLs and CELs. The results are 
summarised in Table 3; the second, third, and last columns give the 
diagnostic lines, their line ratios, and the resulting ri E and T e , re¬ 
spectively. The numbers in the first column indicate the ID number 
of each curve in the 7? e -T' e diagram, presented in Figure 4. The thick 
lines indicate the diagnostic lines for T € , whereas the dashed lines 
correspond to zz E diagnostics. Below, we explain in detail the zz E and 
T € calculations using RLs and CELs. 

3.2.1 RL diagnostics 

The intensity ratio of a high-order Paschen line Pzz (where z? is 
the principal quantum number of the upper level) to a lower-order 



Figure 4. n e -T € diagram based on diagnostic CEL ratios. Each curve is 
labelled with an ID number, listed in Table 3. The thick lines indicate the 
lines of T t diagnostic. The dashed lines indicate the diagnostic lines of zz E . 
The n e -T € curve calculated from the [N i] 7(45198)//(4 5200) is out of the 
plane, because zz E =1010cm -3 . 


Table 3. Plasma diagnostics. 


ID 

n e diagnostic 

Value 

Result 

(cm- 3 ) 

(1) 

[S n](4 6116)1 (A 6731) 

0.60+0.03 

3970+690 

(2) 

[O n](4 3726)/(4 3729) 

1.82+0.13 

3540+790 

(3) 

[S m](/l 18.7 /zm)/(4 33.5 pm) 

2.30+0.19 

6860+1060 

(4) 

[O n](4 3726/29)/(4 7320/30) 

7.42+0.54 

6870+650 

(5) 

[Cl m](/l 5517)/(4 5537) 

0.69+0.03 

8280+870 

(6) 

[Ar rv] (4 47 11 )/(A 4740) 

0.68+0.01 

11320+540 

(7) 

[Ne m](/l 15.6/im)/(/l 36.0jum) 

15.43+1.28 

23 250+8230 


[N i](/15198)/(/15200) 

Paschen decrement 

1.42+0.12 

1010+260 

25 700+6300 

ID 

T e diagnostic 

Value 

Result 

(K) 

(8) 

[S hi] (4 9069)/(,l 18.7/zm+33.5/im) 

0.56+0.05 

11270+1060 

(9) 

[Ne iii](4 3869+4 3967)/(4 15.6/zm) 

1.87+0.12 

12 800+300 

(10) 

[O m](4 4959+4 5007)/(4 4363) 

106.4+1.8 

12 380+80 

(ID 

[N n](4 6548+4 6583)/(4 5755) 

54.4+2.0 

12090+230 

(12) 

[Ar m](4 7135+4 7751)/(4 5191) 

107.8+5.2 

12000+250 

(13) 

[Cl iv](4 8046)/(4 5323) 

26.8+5.7 

11810+1020 

(14) 

[S n](4 6717+4 6731)/(4 4069) 

3.81+0.16 

12 760+810 

(15) 

[S m](4 9069)/(4 6312) 

5.62+0.48 

12 920+620 


He i(4 7281 )/(4 5876) 

0.066+0.003 

13 600+1770 


He i(4 7281)/(4 6678) 

0.145+0.003 

11700+550 


(Paschen Jump)/(P11) 

0.018+0.002 

11690+2090 


Note — Corrected recombination contribution for [O ii] AA 7320/30 A and 
[O iii] ,14363 A lines. 

Paschen line is sensitive to the zz E , as demonstrated in Fang & Liu 
(2011). The Paschen decrement can be an n E indicator for high- 
density regions. As such, we ran small grid models to estimate zr E . 

We calculated the electron temperatures derived from Hei 
lines T E (Hei) using the two Hei line ratios and the emissiv- 
ities of these Hei lines from Benjamin et al. (1999) for the 
case of zz E =10 4 cm -3 . For He +,2+ abundance calculations, we 
adopted the r E (Hei)=12650±1320K, which is the average be¬ 
tween the values derived from Hei 7(47281 A)//(4 5876 A) and 
He 1 1(A 7281 A )/I(A 6678 A). 

We calculated the electron temperature derived from the 


MNRAS 000, 1-19 (2015) 













6 M. Otsuka 


Paschen jump discontinuity T f (PJ) by applying Equation (7) of 
Fang & Liu (2011): 


r e (PJ) = 8.72 x 


He 4 

1 + 0.52-+ 4.40 

H+ 


He 2+ \ 


/.i(pj) r 1 - 77 

/(Pll)j 


( 2 ) 


where He +,2+ /H + is the ionic He + abundance of 9.63(-2) 4 and He 2+ 
abundance of 1.08(-2), respectively, /,i(PJ) is the strength of the 
Paschen jump - i.e. the difference in de-reddened flux density be¬ 
tween 8169A and 8194A in ergs -1 cm -2 A -1 , and /(PI 1) is the 
flux of the Hid8862 A in ergs -1 cm -2 . Because the Hi48862 A 
line is in an order gap, we estimated /(PI 1) using the theo¬ 
retical ratio of /(PI 1 )//(P12)= 1.3 in the Case B 7 E =10 4 K and 
n E =10 4 cm -3 (Storey & Hummer 1995). We adopted the T e (PJ) for 
the RL C 2+,3+ ' 4+ and 0 2+ abundance calculations. 




6820 6825 6830 6835 
Wavelength (A) 


Figure 5. The line profiles of the [Kriv] 4 5867.74 A (left panel) and the 
[Kr iii] 4 6826.85 A (right panel). The grey lines are the observed spectra at 
the rest wavelength, which is shifted by the systemic heliocentric radial ve¬ 
locity (+133.12 km s -1 ). The sky emission lines were removed out as much 
as possible. The red lines are the results of Gaussian fitting. The laboratory 
wavelengths of these lines are indicated by the dashed lines. 


3.2.2 CEL diagnostics 

We derived the CEL T € and n € by solving for level populations 
using a multi-level atomic model (>5), resulting in eight calculated 
n E , including n E ([N i]), as well as eight 7 E . 

For the [O n] 44 7320/30 A and [O m] A 4363 A lines, we sub¬ 
tracted the recombination contamination from the 0 2+ and 0 3+ 
lines using Equations (3) and (4), originally given by Liu et al. 
( 2000 ), 

MO II] AA 7320/30) ( T t f 44 0 2+ 

—m — =9J6 (ios) if’ (3) 

MOIHM43Q) = 124 inr x oy 

/(H/3) \ 10 4 / H+ 

In the 7 r ([Oii] AA 7320/30) calculation, we adopted the r e (PJ) and 
recombination 0 2+ abundances (3.12xl0 -4 , see Section 3.3.2). For 
the l ( ([Om] 44363) calculation, we used the 0 3+ abundance de¬ 
rived from the fine-structure [Oiv] 4 25.89/rm line (6.93xl0 -6 , 
See section 3.3.1) and 7 E ([Nem]). 4([On] 447320/30) and 
/ff(([Om] 44363) were 0.86 and 0.01, respectively. 

Walsh et al. (1997) reported a single 7/ of 12 400±400 K and 
n e of 6000±1500cm - ’. In the present study, we first calculated 
all h e (CEL)s under their derived value / E =12400 K, except for 
n f ([N i]), where we assumed a T € of 10 000 K because the effec¬ 
tive collisional strengths of the Mg° is available at 7/= 10 000 K 
(See Section 3.3.1 for details). The T E =10000K would be a good 
choice for estimates of neutral atoms such as N° and O 0 even 
in high-excitation PNe such as BoBnl (Otsuka et al. 2010). We 
derived 7 E ([Nem]) using n E ([Ariv]) because the ionisation po¬ 
tentials (IPs) of both ions are similar. We derived 7 e ([0iii]) and 
r E ([Cl iv]) under the average n E among n E ([Ar iv],[Cl hi] ,[S iii]); the 
IPs of the [O iii] and [Cliv] are approximately median between 
n E (Cni&[Sm],[Ariv]) or slightly larger than these ions. We cal¬ 
culated 7 E ([Arm]) and two 7 E ([S m]) under the average n € between 
n E ([Clm].[S iii]). 7 e ([Nii]) and 7 e ([Sii]) were calculated under the 
average n E among the two n E ([On]) and n E ([S n]). 


3.3 Ionic abundances 

3.3.1 CEL ionic abundances 

We obtained 21 ionic abundances using CELs, as listed in Table B1 
of Appendix B. We carefully determined the adopted T e and n E 
pairs for the calculations of each ion using the n E -7 E diagram. The 
adopted 7/ and n E values for each ion are listed in the fourth and 
fifth columns of Table Bl, respectively. For the Mg° abundance 
calculation, we adopted the effective collision strengths installed in 
the three-dimensional Monte Carlo photoionisation code Moccasin 
(Ercolano et al. 2003) and the transition probabilities of Mendoza 
(1983). Because the effective collision strengths of Mg° are only 
available at the case of 7 E =10 4 K, we adopted this temperature. 
7 " e = 10 4 K and n E ([Ni]) were used for the N° and Mg° calcula¬ 
tions. The ionic abundances were calculated by solving the sta¬ 
tistical equilibrium equations for more than five levels, given the 
adopted T € and n E , except for Ne + , where we calculated its abun¬ 
dance using the two-energy level model. The last column contains 
the resulting ionic abundances relative to the H + , X m+ /H + , and their 
l-o" errors, which include errors from line intensities, T e , and n € , 
but does not consider the accuracy of the atomic data. In the 0 + 
and 0 2+ abundance calculations, we corrected the line intensities 
of the [O n] 44 7320/30 A and [O m] 4 4363 A, respectively. In the 
last line of each ion’s line series, the adopted ionic abundance and 
its error are indicated in boldface. These values are estimated from 
the weighted mean of the line intensity. 

As presented in Table Bl, the calculated abundances of each 
ion using different transition lines (i.e., fine-structure, nebular, au¬ 
roral, and trans-auroral lines) are well consistent (within the error), 
indicating proper selection of the adopted 7/ and n E pairs for each 
ion and accurate measurement of the flux. 

We report the first detection of the [Kriv] 4 5867.74 A and 
[Kr iii] 4 6826.85 A in Wrayl6-423. The line profiles are presented 
in Figure 5. The sky emission lines between ~6820 and 6830 A 
were subtracted as much as possible. Because there are no can¬ 
didate lines within 0.18-0.20 A from the measured central wave¬ 
lengths of the above lines (5867.74 A and 6827.02 A at rest), we 
concluded that the lines corresponded to krypton. Krypton lines 
have been observed in the Sgr dSph PNe, specifically, BoBnl 
(Otsuka et al. 2010) and Hen2-436 (Otsuka et al. 2011). The Kr 2+ 

4 X(—Y) means X X 10 1 hereafter 
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abundance in Wrayl6-423 is comparable to that in Hen2-436, 
which shows the same metallicity. 

3.3.2 RL ionic abundances 

We list the RL ionic abundances in Table B2 of Appendix B. To 
the best of our knowledge, RL C 3+,4+ and 0 2+ for Wray 16-423 are 
reported for the first time. Because we detected Cn,m,iv and On 
lines, we were able to calculate the elemental C/O ratio using the 
same type of emission lines. 

In the abundance calculations, we adopted the Case B assump¬ 
tion for the lines from levels that have the same spin as the ground 
state, but the Case A assumption for lines of other multiplicities. In 
the last line of each ion’s line series, we present the adopted final 
ionic abundance and its l-cr error estimated from the line intensity- 
weighted mean, indicated by boldface. Because the RLs are n € - 
insensitive under <10 8 cnf 3 , we evaluated the recombination co¬ 
efficients (provided as a polynomial function of T € ) at a constant 
n E =10 4 cm~ 3 for all lines. 

In the He + abundance calculations, using He i A 4471 A and 
He iT5876 A lines, we assumed Case A conditions and adopted 
the recombination coefficients of Pequignot et al. (1991). We sub¬ 
tracted the contribution to these lines by collisional excitation from 
the He 0 2 s 3 S level, using the formulae given by Kingdon & Ferland 
(1995). 

We detected the different multiplet O n lines. Because the elec¬ 
tron density in the 0 2+ emitting region is >10 4 cnT 3 from our 
plasma diagnostics, the upper levels of the transitions in the VI 
O n line would be in local thermal equilibrium (LTE). Therefore, 
we did not correct the VI On line intensities. To calculate the 0 2+ 
abundances, we excluded the VI On,14641.81/4651.33/4673.73 A 
and V2 O ii A 4325.76 A lines. The first two V1 lines are likely con¬ 
taminated with N hi A 4641.85 A and C hi A 4651.47 A, respectively. 
The 0 2+ abundance derived from the On,14673.73 A line is larger 
than that from other VI lines. The last V2 line is likely contami¬ 
nated with C ii A 4325.83 A. 

The RL 0 2+ abundance is larger by a factor of 1.64±0.43 
larger than the 0 2+ CEL. According to Liu (2006), the 0 2+ 
discrepancy factor in Wray 16-423 would be a lower limit. As 
proposed for the halo H4-1 (0 2+ discrepancy factor=1.75±0.36; 
Otsuka & Tajitsu 2013), we speculate that the 0 2+ discrepancy in 
Wrayl6-423 could be explained because the estimated CEL 0 2+ 
abundance could be increased by ~0.2 dex or less if we include 
temperature fluctuations, which is originally proposed by Peimbert 
(1967). Since the 0 2+ discrepancy factor is already low, we do not 
discuss further in the paper. 

3.4 Elemental abundances 

The elemental abundances of the nebula are listed in Table 4. The 
types of emission lines used for the abundance estimations are 
specified in the second column. The fifth column lists the num¬ 
ber densities relative to the solar value, where [X/H] corresponds 
to log! 0 (X/H)-logj 0 (X/H) e . 

To estimate elemental abundances using only the observed 
ionic abundances, we introduced an ionisation correction factor, 
ICF(X), based on the ionisation potential. The adopted ICF(X) for 
each element is listed in Table B3 of Appendix B, and the value is 
listed in the seventh column of Table 4. We calculated the adopted 
ICF(Kr) using Equation (5) of Sterling et al. (2007), based on pho¬ 
toionisation models. 


In the eighth column of Table 4, we list the elemental abun¬ 
dances of Wray 16-423 compiled by Zijlstra et al. (2006a), origi¬ 
nally calculated using the photoionisation model of Dudziak et al. 
(2000) based on the analysis of optical spectra by Walsh et al. 
(1997). The ninth column is the result of Walsh et al. (1997). In 
general, our calculated abundances, except for S and K, are in ex¬ 
cellent agreement with prior results. Walsh et al. (1997) reported 
a C 2+ /0 2+ ratio of 3.4±1.2, derived from Cii and [Om] lines 
(C 2+ (RL)/O 2+ (CEL)=2.48±0.74 in the present work). Assuming 
that C(RL)/0(CEL) = C 2+ /0 2+ , their C abundance would be 5.99(- 
4)±2.17(^1), which corresponds to 8.78±0.16 dex. The differences 
of the S and K abundances between Zijlstra et al. (2006a) and ours 
might be caused by adoption of different T e for these ionic abun¬ 
dances and overestimate of the S 3+ abundance. The [S iv] lines are 
not seen in optical but in mid-infrared wavelength. Their photoion¬ 
isation without dust grains would affect the nebula’s density and 
temperature structure and the fraction of each element’s ion. 

The RL C/O ratio of 2.46±0.80 indicates that Wrayl6-423 is 
presumably a C-rich PN. Elemental RL and CEL O abundances 
seem to follow the trend; as shown in Liu (2006), Wang & Liu 
(2007), and Delgado-Inglada & Rodriguez (2014), the RL C/O ra¬ 
tio seems to be consistent with the CEL C/O ratio within errors 
(.i.e., the RL C/O - the CEL C/O), except for several objects. If 
this is also the case for Wrayl6-423, the CEL C abundance would 
be 8.70±0.14, where the upper value is consistent with the lower 
value of the RL C abundance. 

The noble gas Ar is not easily tied up in dust grains, and is 
also not synthesised in low-intermediate mass stars. Therefore, Ar 
can be used as a metallicity indicator of the progenitor star. The 
[Ar/H] abundance is consistent with that in Hen2-436 (-0.60±0.33, 
Otsuka et al. 2011). Like for Hen2-436, [Fe/H] is significantly 
lower than [Ar/H], Because [Fem] intensities are comparable to 
BoBnl (Otsuka et al. 201 1) and their Fe 2+ abundances in Wrayl6- 
423 and BoBnl are consistent, the small Fe abundance in Wrayl6- 
423 cannot be attributed to inaccuracies in Fe 2+ atomic data or the 
selection of T € for their ionic abundances. Thus, a large fraction of 
Fe atoms would exist as dust grains. 

An important contribution of our work is the detection of the 
light j-process element Kr and calculation of its abundance. This 
required verification of whether or not the Kr abundance value for 
Wrayl6-423 is reasonable. Sterling & Dinerstein (2008) examined 
the relationship between the [Kr/O] and C/O ratio, based on their 
AA-band spectroscopic survey for Galactic PNe, from which the 
following relationship was obtained: 

[Kr/O] = (0.38 ± 0.10) + (0.79 ± 0.29) ■ log 10 (C/O). (5) 

By substituting the observed RL log 10 (C/O) ratio, we obtained 
a predicted [Kr/O] value of +0.69±0.42, which is in good agree¬ 
ment with the observed [Kr/O] of +0.65±0.11 from Kr and O 
CELs. The authors also showed the light ^-enhancement ([ls/Fe], 
Is: Sr, Y, and Zr) versus the C/O ratio among AGB stars and CH 
sub-giants: 

[ls/Fe] = (0.89 ± 0.06) + (1.47 ± 0.18) • log 10 (C/O). (6) 

As described above, the calculated Fe abundance in Wray lb- 
423 may not accurately reflect the correct metallicity. Given that 
[Fe/H] is equal to [Ar/H], the [ls/Ar]=[Kr/Ar] is moderately en¬ 
hanced, +0.82±0.13. This value is slightly smaller than the pre¬ 
dicted [ls/Fe] of 1.46±0.44 from Equation (6), although it depends 
on the adopted metallicity. There are a few results on the .s-process 
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Table 4. Elemental abundances from CELs and RLs. 


(1) 

X 

(2) 

Types of 
Emissions 

(3) 

X/H 

(4) 

log m (X/H) 

+ 12 

(5) 

[X/H] 

(6) 

l°glo(X 0 /H) a 
+ 12 

(7) 

ICF(X) 

(8) 

log 10 (X/H) 

+ 12 (Ref.l) 

(9) 

log 10 (X/H) 

+12 (Ref.2) 

GO) 
log 10 (X/H) 
+ 12 (Ref.3) 

He 

RL 

1.06(—1)±9.56(—3) 

11.02+0.04 

+0.09+0.04 

10.93+0.01 

1.00 

11.03+0.01 

11.03+0.02 

11.09+0.05 

C 

RL 

8.29(—4)±1.57(—4) 

8.92+0.08 

+0.53+0.09 

8.39+0.04 

1.00 

8.86+0.06 


9.08+0.35 

N 

CEL 

4.94(-5)±5.04(-6) 

7.69+0.04 

-0.14+0.07 

7.83+0.05 

24.81+2.12 

7.68+0.05 

7.62+0.10 

7.38+0.22 

O 

CEL 

2.06(-4)±4.10(-6) 

8.31+0.01 

-0.38+0.05 

8.69+0.05 

1.00 

8.33+0.02 

8.31+0.07 

8.35+0.03 


RL 

3.37(—4)±8.90(—5) 

8.53+0.11 

-0.16+0.13 

8.69+0.05 

1.08+0.03 



9.56+0.11 

Ne 

CEL 

4.09(-5)±2.93(-6) 

7.61+0.03 

-0.26+0.10 

7.87+0.10 

1.00 

7.55+0.03 

7.50+0.08 

7.62+0.05 

S 

CEL 

2.51(—6)± 1.58(—7) 

6.40+0.03 

-0.79+0.05 

7.19+0.04 

1.00 

6.67+0.04 

6.48+0.08 

6.83+0.08 

Cl 

CEL 

5.42(—8)±5.66(—9) 

4.74+0.05 

-0.77+0.30 

5.50+0.30 

1.00 

4.89+0.18 


4.86+0.09 

Ar 

CEL 

1.00(-6)±3.85(-8) 

6.00+0.02 

-0.55+0.08 

6.55+0.08 

1.00 

5.95+0.07 

5.88+0.08 

5.95+0.32 

K 

CEL 

2.02(-8)±2.30(-9) 

4.31+0.05 

-0.80+0.07 

5.11+0.05 

3.06+0.23 

4.65+0.22 



Fe 

CEL 

3.62(-7)±5.62(-8) 

5.56+0.07 

-1.91+0.07 

7.47+0.03 

23.81+2.12 



5.73+0.20 

Kr 

CEL 

3.52(-9)±4.88(-10) 

3.55+0.06 

+0.27+0.10 

3.28+0.08 

1.00 



3.62+0.13 


References (1) Zijlstra et al. (2006a), Dudziak et al. (2000), (2) (2) Walsh et al. (1997) for Wray 16-423; (3) Otsuka et al. (2011 ) for Hen2-436. 
“Asplund et al. (2009) for N and Cl, and Lodders (2003) for the other elements. 


in carbon stars in the Sgr dSph. For example, de Laverny et al. 
(2006) reported ^-process enhancements in two C-rich stars 
(C/O=1.05-1.18). The respective [ls/Fe]s in IGI95-C1 ([M/H]=- 
0.8, [ls/M]=+0.6) and IGI95-C3 ([M/H]=-0.5, [ls/M]=+1.0) are 
comparable to that in Wrayl6-423. In comparison with previous 
studies of PNe and AGB stars, we concluded that the Kr abundance 
value obtained for Wrayl6-423 is reasonable. 

In the last column of Table 4, we list the elemental abun¬ 
dance of Hen2-436 by Otsuka et al. (2011); the abundances ob¬ 
tained in the present study, with the exception of N, are close to 
those for Hen2-436. The O and Ar abundances indicate that both 
PN progenitors would form under a similar evolutionary stage of 
the Sgr dSph. Dudziak et al. (2000) pointed out that the nearly iden¬ 
tical abundance patterns for these two Sgr dSph PNe provide evi¬ 
dence that the progenitors formed in a single starburst event (~5 
Gyrs ago, Zijlstra et al. 2006a) within a well mixed ISM. The age- 
metallicity relationship for the Sgr dSph globular cluster candidates 
by Law & Majewski (2010) indicates that the [Fe/H]—0.6 corre¬ 
sponds to ~8 Gyrs ago. Layden & Sarajedini (2000) reported that 
the Sgr dSph field is composed of several different populations; 
the [Fe/H]=-0.7±0.2 population is 5±1 Gyrs old and the [Fe/H]=- 
0.4±0.3 population is 0.5-3 Gyrs old. Although the exact [Fe/H] 
abundances in Wray 16-423 and Hen2-436 are unknown, the pro¬ 
genitor stars of Wray 16-423 and Hen2-436 most likely formed 0.5- 
6 Gyrs ago. 


3.5 [WC] type central star 

In Figure 6, we present the profile of the C iw AA 5801/12 A 
band composed of narrow nebular Civ AA 5801/12 A lines and a 
weak broad wing component (perhaps originating from a stel¬ 
lar wind). We fitted the band profile by three Gaussian compo¬ 
nents: in Figure 6, the thin black lines are the components of 
the nebular C iv AA 5801/12 A, the red line indicates our interest¬ 
ing broad component, and the thick black line is the sum of these 
three. The FWHMs of the Civ AA 5801/5812 A lines are 0.75±0.04 
and 0.83±0.05 A and that of the broad component (centered at 
5805.96 A in the rest frame) is 32.12±2.15 A, where its FWHM ex¬ 
pansion velocity is 1660± 110 km s -1 . The CivAA 5801/12 A band 
profile in Wrayl6-423 is similar to the Sgr dSph PN StWr2-21; 
Kniazev et al. (2008) de-convolved the C iv AA 5801/12 A into three 
components, which included a broad component having a FWHM 
of 42.3±3 A. Although the nebular C iv AA 5801/5812 A lines seem 



Wavelength (A) 


Figure 6. The observed profile of the Civ AA 5801/12 A band indicated by 
the grey line. The observed wavelength is corrected by the heliocentric ra¬ 
dial velocity of + 133.12 km s . An interesting feature (perhaps originating 
from a stellar wind), is indicated by the red line. 


to be completely embedded in the broad component in Hen2-436, 
the FWHM of the C rv AA 5801/12 A band is 32.8±0.2 A, measured 
from the same ESO VLT/FEROS2 spectra used in Otsuka et al. 
(2011), and consistent with the findings of Walsh et al. (1997, 
31.1 A). We did not detect any other broad components, such as 
O vi AA 3811/34 A. 

A consensus has not been reached regarding the classifi¬ 
cation of the CSPN of Wrayl6-423. Walsh et al. (1997) classi¬ 
fied it as a [WC8] or a weak emission line star (WELS), on the 
basis of the strengths of the Cni/14650A (FWHM=9.9A) and 
CIV/M5801/12A (FWHM=14.7A) bands, using low-dispersion 
spectra (R could be ~1200 around 5800 A); the broad component 
under the nebular Civd/15801/5812 A lines was not detected. The 
broad component of the Cm and Crv complex centered around 
4650 A, originating from the stellar winds, is not clearly seen in 
our HDS and FEROS spectra. There is no broad component in the 
C iiT4267 A, too. 

Taking into account that their estimated CSPN temperature 
(85 000 K) is hot for a [WC8] PNe, Walsh et al. (1997) concluded 
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Wavelength (fim) 


Figure 7. Spitzer/IRS spectrum of SMP LMC19 in the LMC. The intensity- 
scaled spectrum of Wrayl6-423 is indicated by the red line. The scaling 
factor is 2.98, for adjustment to the 16-36//m flux density of LMC19. 


that Wrayl6-423 could fall into WELS. Acker et al. (2002) anal¬ 
ysed 42 emission line nuclei of PNe, and reported that the [WC4] 
CSPNe show the highest temperatures among the [WC] class 
(54950-91 200 K). The determination of CSPN temperatures can 
be difficult, as pointed out by Acker et al. (2002), which may ex¬ 
plain why the CSPN temperatures of the [WC] class PNe have a 
relatively wide range. Even after the CSPN temperature was cor¬ 
rected up to 107 000 K by photoionisation models based on the low- 
dispersion spectra of Walsh et al. (1997), Gesicki & Zijlstra (2003) 
and Zijlstra et al. (2006a) supported the conclusion by Walsh et al. 
(1997). 

According to the classification of Acker et al. (2002), the 
broad FWHM of the stellar C rv component suggests that the CSPN 
of Wray 16-423 may be a [WC4]-type central star. 


3.6 PAHs and dust features 

Wray 16-423 shows PAH bands and dust features, specifically, 
6-9 pm and 10-14yum PAH bands and broad \\ pm, 16-24 pm, 
and 30yum features. These PAH bands and these dust fea¬ 
tures have been found in C-rich MC PNe (Stanghellini et al. 
2007; Bernard-Salas et al. 2009; Garcfa-Hernandez et al. 2012; 
Matsuura et al. 2014; Sloan et al. 2014). Bernard-Salas et al. 
(2009) reported broad 11 gm, 16-24/rm 5 , and 30 gm features in 8 
C-rich MC PNe. Because the typical LMC metallicity is close to 
that of Wray 16-423, LMC PNe provide a good comparison in terms 
of dust production in a metal-deficient environment. Therefore, we 
looked at the LMC PN sample to compare the mid-IR spectrum of 
Wray 16-423 with those of selected LMC PNe. 

Accordingly, we found that SMP LMC19, 25, 48, and 78 ex¬ 
hibit broad 11/16-24/30 gm features, comparable to the dust feature 
strengths measured in Wray 16-423. Among them, the mid-IR spec¬ 
trum ofLMC19is most similar to that of Wray 16-423 (Fig. 7); the 
strengths of the broad 16-24 gm and 30 pm features and the global 
~5-37 gm SED are similar as well. 

For a more detailed look at these PAH and broad dust fea¬ 
tures, we subtracted the expected local dust continuum, determined 
by fifth-order spline-function fitting. The resultant spectra are dis- 


5 Bernard-Salas et al. (2009) call this Bump 16. 



5 10 15 20 25 30 35 

Wavelength (yu.m) 


Figure 8. Spitzer/IRS spectra of SMP LMC25,48, 19, and 78 and Wray 16- 
423 (Wrayl6). The local continuum was subtracted by fifth-order spline- 
function fitting. Flux densities are normalised and shifted to emphasise the 
16-24 and 30//m features. 


played in Figure 8. The flux densities are scaled and shifted to em¬ 
phasise broad 16-24yum and 30 pm features. 

In the following sections, we describe the 6-9 pm and 10- 
14/rm PAH bands, the broad llyum, 16-24/jm, and 30 gm band 
features in Wrayl6-423, and LMC19, 25, 48, and 78. Dust and gas 
mass estimates for these LMC PNe will be discussed in a future 
study. 


3.6.1 The 6-9pm PAH band 

The classification of the 6-9 pm PAH band is based on the positions 
of 6.2, 7.7, and 8.6/jm features at the intensity peaks. According to 
Peeters et al. (2002) and Matsuura et al. (2014), the 6-9 gm PAH 
profiles have peak positions near 6.22-6.3/rm, 7.6-8.22gm, and 
>8.6 pm. Our classifications are summarised in Table Cl of Ap¬ 
pendix C. From the overview of Peeters et al. (2002), both the 6.2 
and 7.7 pm PAHs in Wrayl6-423 are Class B, with corresponding 
peaks at 6.24±0.01 pm, and 7.78±0.04yum, respectively. The 6.2 
and 7.7 pm PAHs in our comparison LMC PNe fall into Classes B 
and /ft, respectively. The 8.6/rm PAH in Wray 16-423 could not be 
resolved, due to a low signal. The S.6pm PAHs in LMC PNe are 
classified as Class B. 
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Figure 9. The 9.5-13.5fim emission complex showing the 10-14pm PAH 
bands and the broad 11 pm feature. The vertical dashed lines mark the 
positions of (a) [Siv] /110.51/im, (b) 11.3 jam PAH, (c) 12pm PAH, (d) 
Hid 12.38pm. and (e) the 12.7 pm PAH + [Nen] 4 12.80pm complex. The 
grey line shows an 750/SWS spectrum of the Galactic C-rich AGB star W 
Ori, with the 11 pm SiC feature. 


3.6.2 The 10-14 /.un PAH band and the broad 11 [im feature 

Figure 9 shows the 10-14 pm PAH bands complex and the broad 
11 pm feature. The flux densities are normalised with respect to the 
11.3-pm PAH emission peak. The positions of the PAHs and the 
atomic gas emission lines are shown as vertical dashed lines with 
the lowercase letters a-e. 

In the study on PAHs in LMC post-AGB stars, Matsuura et al. 
(2014) classified the 10-14pm PAH feature into four types: a, /3 , 
y, and 6, depending on the peak wavelengths and the shapes of 
the sub-features. The Class a (including the broad 11 pm feature) 
shows a sharp feature at 11.3 pm, a well-isolated feature at 12.7 pm, 
and a very weak feature at 12pm. The 11.3pm PAH emission 
in Class a is the strongest among the three features. The respec¬ 
tive features would be from the solo (11.3pm PAH), duo (12 pm 
PAH), and trio (12.7 pm PAH) out-of-plane C-H bending mode 
of PAH (Sloan et al. 2014). We classify the 10-14pm features in 
Wrayl6-423, LMC19, and LMC78 as Class a. The flux density 
for these three PNe reaches its maximum in this range at the peak 
of the 11.3 pm PAH emission (central wavelengths: 11.28, 11.30, 
and 11.29pm, respectively); 12pm PAHs are also observed in all 
three. The 12.7 pm PAHs are observed in LMC19 and 78, while 
the presence of this PAH in Wray 16-423 seems to be merged with 
[Nen] A 12.80 pm. 

LMC25 and 48 are characterised as Class <5, by the broad fea¬ 
ture extending from 10 to 14 pm with the 11.3 pm PAH on the 



Wavelength (pm) 


Figure 10. The 14-37.5 pm continuum subtracted spectra, showing the 
broad 16-24pm and 30pm features. The vertical dashed lines mark 
the positions of (a) [Ne iii] 4 15.56pm, (b) and (c) 16.4pm PAH and 
15-18pm plateau (possible), (d) [S m] ,118.71 pm, (e) [Ne iv] ,124.32pm, 
(f) He lid 25.31 pm, (g) [O iv] 4 25.89pm, (h) [S iii] 433.50pm, and (i) 
[Ne iii] A 36.01 pm. 


top (Matsuura et al. 2014). Class 8 spectra are quite different from 
Class a spectra, in terms of the underlying broad 11 pm feature and 
its central position. We fitted Gaussians to the PAH and atomic gas 
emission lines, and subtracted the fits from the spectra. We then fit¬ 
ted the broad 11 pm feature in the residual spectra was fitted using a 
single Gaussian component to determine their central wavelengths 
A c and FWHMs (Table Cl ). The average peak position and FWHM 
of the broad feature are 11.92±0.05pm and 1.74±0.11 pm, respec¬ 
tively, for Class a, while those in Class 6 are 11.46±0.01 pm and 
1.86±0.03 pm, respectively. 

Matsuura et al. (2014) argued that the broad 11pm feature 
characteristic to Class S should be ascribed to silicon carbide (SiC). 
They used fits for the 10-14 pm spectra of three post-AGB stars 
using the NASA/Ames PAH IR-spectral database (Boersma et al. 
2014; Bauschlicher et al. 2010) 6 and the /SO/SWS spectrum of the 
C-rich AGB star W Ori as the SiC template (see also Fig. 9). 
Sloan et al. (2014) analyzed the broad 11 pm features and associ¬ 
ated PAHs in LMC and SMC post AGB stars. They call the band 
profiles similar to LMC25 and LMC48 “Big-11”. They demon¬ 
strated that this Big-11 is composed of SiC and PAH, and they con¬ 
cluded that the Big-11 feature is due primarily to SiC (88 % of the 
total flux). Otsuka et al. (2015) argued that the broad 11 pm fea¬ 
ture in W Ori could be reproduced using the absorption efficiency 
Q a of a spherical a-SiC grain. Assuming that SiC grains greatly 
contribute to the 10-14 pm feature of Class 6 , the wider FWHM in 
Class 8 compared with Class a would then be due to the contribu¬ 
tion from the 12 and 12.7 pm PAHs. 
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Figure 11. The plot between the 6-9 /im PAH and the broad 16-24//m fea¬ 
ture strengths. 


3.6.3 The broad 16-24 fim feature 

Broad 16-24/mi and 30/mi features are presented in Fig. 10. 
Bernard-Salas et al. (2009) and Garcla-Hernandez et al. (2012) re¬ 
ported the detection of this feature in MC PNe. Matsuura et al. 
(2014) detected this feature in the LMC post-AGB star IRAS0588- 
6944. In our galaxy, a broad 16-24/im feature has been found 
in PNe Ml-11 (Otsuka et al. 2013), K3-17, M2-43 (Hony et al. 
2002b), and the proto-PN IRAS01005+7910 (Zhang et al. 2010). 

The 16-24 gm feature is narrower in LMC25, 48, 78 (See Ta¬ 
ble Cl); this could be due to the contribution of the blue wing of 
the 30 gm feature. The band profiles in LMC25, 48, and 78 are 
slightly asymmetric; the /l c s of the 16-24/mi feature in these three 
are blue-shifted, compared with those of Wrayl6-423 and LMC19. 
The asymmetric profile could be due to a small plateau superposed 
on the broad 16-24/jm around 15-18/mi. The 15-18/mi plateau 
seems to be possibly in LMC19 (/l c ~17.10/mi), but its presence is 
unclear in Wrayl6-423. In Figure Cl of Appendix C, we compare 
our sample PNe with the Galactic PN K3-60 (AORKEY: 19903744, 
PI: H. Dinnerstein, indicated by the grey lines). We use the K3- 
60 spectrum to explain the position and width of the 15-18/mi 
plateau, as an example. K3-60 shows the 15-18/mi plateau but no 
broad 16-24/ini feature. The 6-9 and 10-14/rm PAH band profiles 
in K3-60 are well matched to those in LMC78, 19, and Wrayl6- 
423. The complex of the 15-18/mi plateau and the 16.4/mi PAH 
also contribute to the broad 16-24/mi feature; however, its con¬ 
tribution seems to be negligibly small in our sample PNe. Even 
if the 15-18/mi plateau plus the 16.4/mi PAH are the major con¬ 
tributors to the 16-24/mi feature, this complex cannot explain the 
asymmetry and difference in width between the observed 16-24 /mi 
features. The red and blue wings of the 16-24 gm feature have dif¬ 
ferent slopes between LMC25/48/78 and LMC 19/Wray 16-423, and 
this cannot be explained by the presence of the plateau only. Note 
that also the broad 30 /mi feature is clearly broader and redder in 
LMC19 and Wrayl6-423. The different widths and positions might 
reflect the colder dust and/or larger dust grains. 

The carrier of the broad 16-24 /mi feature has yet 


6 http://www.astrochem.org/pahdb/ 


Table 5. Fitting results for the broad 30//m feature. The uncertainty of 
7)/(max) is within 0.6 K. 


Fitting 

Fitting wavelength 
(//m) 

T d ( max) 
(K) 

Fitl 

14-22.65 

145.5 

Fit2 

14-65 

131.3 

Fit3 

27-65 

108.1 

Fit4 

14-15, 27-65 

133.0 


to be identified. Van Kerckhoven et al. (2000) noted that the 
6.2/7.7/8.6/11.3/tm PAH emission is from relatively small PAH 
molecules (~50 C-atoms). Additionally, they argued that the 16- 
20/mi plateau - it should be attributed to the 15-18/mi plateau as 
described above - is due to the destruction of large PAHs and PAH 
clusters, containing ~2000 C-atoms. 

Can the 16-24/mi feature be also attributed to large/cluster 
PAHs and the 6-9 /mi PAH form by the destruction of large PAHs 
and PAH clusters? If so, then an anti-correlation may exist between 
the intensities of the 16-24/mi feature and the 6-9 /mi PAH. In 
7 LMC PNe listed in Table C2 of Appendix C and Wray 16-423 
showing both features, we measured the integrated fluxes of the 
6-9/mi PAH band (F(6-9/mi PAH)), the 16-24 /mi feature (F( 16- 
24/mi)) in the local continuum subtracted spectra, and the inte¬ 
grated flux between 5.75-36.8/mi (F(IR)). We define the ratios of 
the F(6-9 /mi PAH)/F(IR) and the F( 16-24 /mi PAH)/F(IR) as fea¬ 
ture strengths, and we present between these strengths in Fig. 11. 
The correlation factor is -0.50, depending on the data of Wray lb- 
423 because the factor is down to -0.16 without Wrayl6-423. 
Therefore, at this moment we cannot conclude that there is a corre¬ 
lation between the 16-24/mi feature and the 6-9/mi PAH band (or 
PAHs themselves). 

Hony et al. (2002b) proposed iron sulfide (FeS) to explain the 
16-24/rm feature. We note Zhukovska & Gail (2008); ? FeS dis¬ 
plays a broad feature centered at 23 /im and two sharp resonances 
at 34 and 38/mi. In our sample PNe, the T c s of the 16-24/mi fea¬ 
ture are not in accordance with that of the broad component of FeS 
at the central wavelength (~23 /mi), and none showed the 34 /mi 
resonance. Thus, Fe depletion in Wray 16-423 does not reflect the 
presence of FeS as a carrier of the 16-24/im feature; however, it 
may imply those of other Fe-based grains. 

Amorphous silicates have broad emission features at two 
peaks near 9.7/mi (Si-0 stretching mode) and ~18/mi (O-Si-O 
bending mode). Thus, the fossil silicate dust that forms in early 
AGB outflow is of interest. Dual-dust PNe have been found in 
the Galactic bulge (e.g., Gorny et al. 2010; Perea-Calderon et al. 
2009). Gorny et al. (2010) reported four PNe showing amorphous 
silicate, PAHs, and crystalline silicate features at 23.5, 27.5, and 
33.8/mi, respectively. Crystalline silicates are not detected in our 
sample. As Figure 7 of Bernard-Salas et al. (2009), the central 
wavelength of the 18/mi bump is not in accordance with that of the 
16-24/mi feature. If amorphous silicate is responsible for the 16- 
24 /mi feature, then a broad bump should be observed in the blue 
shoulder of the broad 11 /mi feature. This is not the case. Therefore, 
the broad 16-24/mi features in our sample cannot be attributed to 
amorphous silicate. 


3.6.4 The broad 30 fim feature 

Our PNe display a broad 30/mi feature, which has been seen 
in many C-rich objects. Hony et al. (2002a) performed a com¬ 
prehensive analysis of the 30/mi feature showing objects with 
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Figure 12. Reproduced SEDs by the Equation (7) (black, green, red, and 
blue lines) overlaid on the Spitzer /IRS spectrum of Wray 16-423 (grey line). 
The point indicated by the diamond shows the upper-limit flux density at 
AKARI/FIS 65 pm. Fits 1-4 are different in terms of the fitting wavelength 
range. The fitting results are summarised in Table 5. 


ISO/ SWS spectra. Due to the wavelength coverage, our measured 
FWHMs were narrower than the measurements of Hony et al. 
(2002a, ~10/jm). Hony et al. (2002a) argued that the feature origi¬ 
nated from a CDE MgS. The CDE MgS has been the major candi¬ 
date for the 30 pm feature. However, as we introduced earlier, the 
origin of this feature is unclear. 

Messenger et al. (2013) performed a detailed analysis on the 
30/rm and the 11 pm SiC features; the authors concluded that the 
30 pm feature is likely due to a carbide or sulfide component, and 
that the abundance of the carrier of the 30 pm feature is linked to 
SiC abundance. This conclusion seems to agree with Sloan et al. 
(2014), who argued that the MgS feature strength climbs just as the 
SiC strength drops, as the dust shells around the stars become red, 
based on the assumption that the SiC/amorphous carbon is coated 
by MgS as the dust grains move away from the central star. Labo¬ 
ratory experiments by Lodders & Fegley (1995) revealed the con¬ 
densation temperatures of different C/O ratios and pressures; for 
the case of C/O = 2.0 and a total pressure of 10~’ bar, the conden¬ 
sation temperatures of graphite, SiC, and MgS grains are 1978 K, 
1736 K, and 1152 K, respectively. 

We measured the F(11 pm SiC)/F(IR) for LMC PNe LMC 8 , 
25, 48, and 85 (Table C2) showing SiC. Among these PNe, an anti¬ 
correlation relationship exists between the SiC and 30 pm intensi¬ 
ties (correlation factor=-0.88), supporting the idea of Sloan et al. 
(2014). We investigated the effective temperatures of 14 LMC 
PNe using Cloudy. The average effective temperature of the four 
SiC-containing PNe (LMC 8 , 25, 48, and 85) is 50750 K, while 
a value of 123 000 K is determined among the non-SiC contain¬ 
ing PNe (LMC 19, 78, and 99, and Wrayl6-423). The CSPNe of 
SiC-containing PNe seem to be cooler and younger than those of 
the non-SiC-containing PNe. Thus, the MgS grains in non-SiC- 
containing PNe may be far from the central stars. In comparison, 
the average A c s of the 30/rm feature in SiC-containing four PNe is 
30.50 pm and that in the remaining non-SiC-containing four PNe is 
31.73 pm. 

However, Zhang et al. (2009) cast doubt on this identification 
for the 30yum feature. They demonstrated that the MgS mass of 
(0.96-7.16)x 10 ~ 5 M 0 , formed from the available S atoms, could 
not account for the observed strength of the 30yum feature in the 


proto-PN HD56126, even if all the dust grains existed as MgS (the 
required MgS mass is >3.5xl0~ 4 M o ). 

CDE iron-rich sulphide, such as Mgo. 5 Feo. 5 S (Fe50S), also 
show a broad feature around 30yum, similar to the CDE MgS fea¬ 
ture (Messenger et al. 2013). In PNe, Henry et al. (2012) reported 
no correlation between the S abundance deficiency and the occur¬ 
rence of the broad 30y^m feature, implying that iron-rich sulfide 
rather than MgS may be favored as a carrier of the 30 pm feature. 
The type of Mg-Fe-S grains that form depends on the amount of 
Mg, Fe, and S atoms in the solid phase. 

Another possible candidate for the 30/rm feature is graphite, 
as proposed by Jiang et al. (2014). The opacity of graphite by 
Draine & Lee (1984) exhibits broad emission (FWHM>60/rm) 
peaked at ~ 30-40 pm. The theoretical opacity of graphite with a 
d.c. conductivity of 100 fi _1 cm 1 by Jiang et al. (2014) shows a 
well-consistent FWHM of ~10/rm to the observed FWHMs re¬ 
ported by Hony et al. (2002a). Otsuka et al. (2014) found that the 
relative strength of the 30yum feature to the underlying continuum 
in fullerene Cgo-containing PNe is constant; they also succeeded 
in fitting both the 30pm feature and the 13-160/rm continuum 
using the synthesised absorption efficiency <2abs(4), derived from 
the combined data of ISO/SWS/LWS and Spitzer/IRS spectra of 
C 6 o-containing PN IC418. Their derived opacity curve is similar 
to that of graphite by Draine & Lee (1984), with the exception of 
the feature strength and FWHM. Otsuka et al. (2014) concluded 
that the 30 pm feature and featureless continuum may be related 
to hydrogenated amorphous carbon (HAC) or graphite. HAC has 
two resonances centered at ~20 pm and ~30 as demonstrated in 
Fig. 9 of Otsuka et al. (2013). However, since its optical constant 
is strongly variable for different chemical and physical conditions 
such as hydrogen-content and grain size, one might be able to re¬ 
produce the strong 30 pm feature plus the weaker sub structures. 

We utilised the method of Otsuka et al. (2014) to verify our 
ability to fit the >14pm SED of Wrayl6-423 with their <2abs(T) 
even if it will not be possible to fit the broad 16-24 pm feature in 
this way. The model of Otsuka et al. (2014) assumes that the dust 
density, as a function of the distance from the CSPN r, is distributed 
around the CSPN with a power-law (oc r~ p ) and that the dust tem¬ 
perature distribution T d (r) also follows a power-law (oc r~ q ). The 
thermal radiation from spherical grains can be estimated by the fol¬ 
lowing: 

J ~T d { max) 

T d a B,(T d ) dT d , (7) 

D(min) 

where C is a constant value, T d ( max) and T d {mm) are the maxi¬ 
mum and minimum temperatures of the 30pm feature, respectively, 
cr=(3 - p)/q > 0, and B d (T d ) is the Plank function. We chose both 
p and q to be 2. The fit is not sensitive to T/ffmin); thus, we let 
7//(min)=20 K, a temperature more or less characteristic of dust in 
the ISM. 

We fit the Spitzer/IRS spectrum for wavelengths longer than 
\4pm by the Equation (7). As a constraint to set the upper limit 
flux density over 36.5 pm by the model fit, we used an upper limit 
flux density / y (65/rm)=88.3 mjy measured from the AKARI/FIS 
Far-infrared All-Sky Survey Maps. We performed four fits, as sum¬ 
marised in Table 5. The fitting regions of each differed. The re¬ 
sultant SEDs are plotted on the Spitzer/IRS spectrum in Figure 12. 
Fitl - a test whether the 14-22 pm dust continuum is linked with the 
30 pm feature - predicts ~ 1.2 times larger flux density at the peak of 
the 30 pm feature. Fit2 converging at 14-65 pm gives a better result; 
however, it cannot explain ~20-24yum of the broad 16-24 pm fea- 
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Figure 13. Radial profile of the hydrogen density adopted in SED modeling. 


Table 6. Best fitting parameters and the derived properties by the Cloudy 
model. 


Parameters of the CSPN 

Values 

L, 

4753 L s 

Teff 

110 110 K 

M v 

2.691 

log g 

6.0 cm s -2 

Distance 

24.8 kpc 

Parameters of the Nebula 

Values 

Abundances 

He: 11.02, C:8.86, N:7.69, 0:8.22, Ne:7.50, 

(l°g 10 (X/H)+12) 

S:6.40, 0:4.71, Ar:5.92, K:4.31, Fe:5.57 
[Si,Mg/H]=-2.00, Others:[X/H] =-0.65 

Geometry 

Spherical 

Shell size 

R m =0.20" (0.02 pc), R out =1.3" (0.16 pc) 

n H 

See Fig. 13 

filling factor 

0.85 

lOgRHyS) 

-11.705 erg s _1 cm -2 (de-reddened) 

m g 

0.59 M 0 

Parameters of the Dust 

Values 

Composition 

3.1 % PAHs and 96.9 % graphite (mass frac.) 

Grain size 

0.0004-0.01 \jim for PAH 

0.005-0.25 lira for graphite 

T d 

113-305 K (PAHs), 43-194 K (graphite) 

md 

1.32(-4)M 0 in total 

m!m g 

2.26(-4) 


ture or the 30/jm. Fit3, focusing only the 30 pm feature, provides 
a good fitting for the 30 gm feature, although it under-predicts the 
flux density below 27 /tm. Fit4 skips the broad 16-24/rm feature 
and gives results similar to those obtained in Fit2. 

Our simple analysis demonstrates that the broad 16-24/tm in 
Wray 16-423 is unrelated to the 30/rm feature. If we focus on the 
27 yum or longer wavelength data, we can reproduce part of the ob¬ 
served 30/rm feature with the synthesised Q d b s (/l) of Otsuka et al. 
(2014). More fits to the 30/rm feature showing PNe (but without 
broad 16-24 /im feature) are necessary to synthesize the absorption 
efficiency of the 30 gm feature for astronomical objects and to sepa¬ 
rate the 16-24 gm feature from the complex of 16-24 gm and 30 yum 
features. 


3.7 SED modeling 

We constructed an SED model to investigate the physical condi¬ 
tions of the gas and dust grains, in order to derive their masses 
using Cloudy clO.OO. 


3.7.1 Model approach 

We set the distance to Wrayl6-423 to 24.8 kpc for comparison 
with the observed fluxes and flux densities. For the incident SED 
from the CSPN, we adopted the model grids of non-LTE line- 
blanketed, plane-parallel, hydrostatic model atmospheres of halo 
stars ([Z/H]=—1) by Rauch (2003). Because the Hen line-fluxes 
are sensitive to the surface gravity logg, we tested small gas- 
only model grids with varying logg that matched the flux of the 
Hen44686A. In the small grids, we set the CSPN’s effective 
temperature to T eff =110 660K and luminosity L, to 4683 L 0 us¬ 
ing Equations (3.1) and (3.2) of Dopita & Meatheringham (1991), 
which are established amongst Magellanic Cloud PNe, assuming 
optically thick nebula. Accordingly, we set logg to 6.0 cm s -2 . 
We kept this surface gravity in the gas+dust full SED model, and 
adopted the above T c t r and L, as a first guess. 

For gas-phase elemental abundances, we adopted the observed 
values listed in Table 4 as the first guess, and refined these to match 
the observed line intensities of each element. We considered the 
CEL O line fluxes to set the O abundance. We adopted the same 
transition probabilities and effective collision strengths of CELs 
used in our plasma diagnostics and abundance determinations. Be¬ 
cause Cloudy does not include any Kr lines, we did not consider 
this element. The abundances for non-calculated elements, except 
for Mg and Si, were held constant at [X/H]=-0.65, which is the 
average between the observed [Cl/H] and [Ar/H], Because the de¬ 
ficiency of refractory elements Mg and Si is unknown, we adopted 
[Mg,Si/H]=-2. 

We determined the radial H-density profile of the nebula, 
based on the radial intensity profile of the HST/WFPC2 F656N im¬ 
age using an Abel transformation under the assumption of spherical 
symmetry. The ionization-bound models over-predicted the [N n] 
and [On] lines and under-predicted the [Om] lines, whereas the 
matter-bound models gave better fit to the observed fluxes of these 
emission lines. Therefore, we assumed the matter-bound condition. 
The outer radius R out was set to 1.3" (0.16 pc) and the inner ra¬ 
dius R in was 0.2" (0.023 pc), from the HST image. We adopted a 
constant filling factor of 0.85. The adopted radial profile of the hy¬ 
drogen density is presented in Figure 13. Our set hydrogen density 
profile appeared to be similar to the photoionisation models without 
dust grains by Dudziak et al. (2000), who adopted the constant hy¬ 
drogen densities of 9500 cnr 3 for 0.023-0.051 pc (covering factor 
=0.17) and 3600 cm- 3 for 0.023-0.08 pc (0.83). 

We attempted to fit the observed SED in the range from 
0.1516 pm {GALEX- Fuv) to 65 pm (AKARI/FIS), assuming that the 
infrared excess around 5-40 yum was due to the emission from PAH 
molecules and carbon-based grains. There are no suitable optical 
constants for the broad 16-24yum and 30-yum features; thus, we 
skipped these bands. The 16-24 pm and 27.4-36.5 yum SEDs were 
out of the fitting range in the model. Here, we adopted graphite 
and PAH. We assumed a spherical shape for both the graphite and 
PAH. In graphite, the spheres were assumed to be randomly ori¬ 
ented, and the ’T/3-2/3” approximation was utilised (for validity of 
this approximation see, Draine & Malhotra 1993). The optical con¬ 
stants were taken from Draine & Li (2007) (PAH-Carbonaceous 
grains) for PAH and Draine & Lee (1984) for graphite. For PAH, 
we adopted the radius a in the range from 0.0004 (30 C-atoms) 
to 0.0011 yum (500 C-atoms) with the standard interstellar dust- 
grain size distribution by Mathis et al. (1977) (i.e., a' 3 5 ). For the 
graphite, we adopted the same size distribution with a=0.005- 
0.25 gm. 

To verify the model fitting accuracy, we evaluated the chi- 
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square value (x 2 ) from the 44 gas emission fluxes, 9 gas-phase 
abundances, 7 broad band fluxes, and the 9 de-reddened flux den¬ 
sities I v from GALEX, W.S7/WFPC2, 2MASS, and WISE bands. 
We used the AKARI/FIS / v (65 pm) to set the dust continuum above 
36.5/rm. 

3.7.2 Comments on SED model results 

The properties of the CSPN, nebula, and dust/molecules matching 
the observed quantities are listed in Table 6. The predicted fluxes 
relative to the H/3, the flux densities at the interest bands, and el¬ 
emental abundances compared with observations are listed in Ta¬ 
ble 7. The reduced chi-square value indicates that our SED model 
can reproduce the observed quantities. 

We may have overestimated the [Ne n] A 12.80/im in the low- 
dispersion Spitzer/IRS spectrum, because this line would be con¬ 
taminated by the 12.7 pm PAH (see Fig. 9). Although we could not 
reproduce these line fluxes within estimated errors, our model is 
globally successful in explaining the observations. Although with¬ 
out the UV spectra we cannot say for certain, the large difference in 
GALEX- Fuv and Nuv might be due to the adopted extinction value 
and reddening correction function for GALEX wavelength. 

The model predicted M v of the CSPN (2.685) is con¬ 
sistent with the observed absolute magnitude of 2.770±0.2 in 
the case of 24.8 kpc within error. The uncertainty of T e u is 
2 000K. Taking into account the distance uncertainty of 0.8 kpc 
(Kunder & Chaboyer 2009), our derived L, of 4753±1230L o 
is fairly consistent with Dudziak et al. (2000), who reported 
4350± 150 L 0 and r eff = 107 000± 10 000 K. The uncertainties of the 
predicted elemental abundances are within 0.1 dex. 

Taking the absolute Spitzer/IRS flux calibration uncer¬ 
tainty -17% (Decin et al. 2004), the dust mass m d would be 
(1.32±0.23)xl0~ 4 Mq. The uncertainty of the dust mass fraction 
would be within ±1,%. The estimated dust mass here would be 
the lower limit value, because we skipped the broad 16-24/jm and 
30/rm bands. If one of the carriers for the 30 pm broad band is the 
Fe-rich sulfide, then the dust mass and the dust-to-gas mass ratio 
could increase. Indeed, the large Fe depletion with respect to Ar is 
confirmed; thus, the missing Fe atoms would exist as Fe dust grains. 
In the next section, we discuss how much dust mass can be added 
if the broad 30 pm feature is from Fe-rich sulfide or MgS. The gas 
mass m g (the sum of the ionised and atomic gas mass within the 
volume occupied within R out ) is 0.59±0.10M Q . We also ran the 
model for the case of 7? out =1.0" (0.12 pc), and we obtained the 
md=(1.05±0.18)xl0~ 4 M 0 and m s =0.45±0.08 M 0 . The estimated 
m g for our R olIt =1.0" and 1.5" models is consistent with the result 
of 0.4 M 0 by Zijlstra et al. (2006a). The ionised gas mass exceeds 
the prior estimate of 0.248 M 0 by Dudziak et al. (2000), which we 
attribute to the difference in H/3 flux and outer radius (log /(H/3)=- 
11.89 erg s _1 cm -2 and R ollt =0.08 pc at 25 kpc). If we adopted their 
/(H/3). then m g is ~0.29 Mo in the case of 7? out =0.12 pc. 

The dust mass in Wrayl6-423 is similar to that in Hen2-436 
((2.9-4.0)xl0~ 4 M 0 , Otsuka et al. 2011), whereas the gas mass is 8- 
12 times larger. Otsuka et al. (2011) and Dudziak et al. (2000) es¬ 
timated 0.05-0.07 M 0 (ionised + neutral gas) and 0.04 Mq (ionised 
gas) for Hen2-436, respectively. Dudziak et al. (2000) argued that 
the outer boundary in Hen2-436 is determined by radiation (ionisa¬ 
tion bound). Indeed, Hen2-436 is very compact; its outer radius 
is -0.03 pc. The small 0 2 7(0 + +0 2+ ) ratio in Hen2-436 (0.85, 
Otsuka et al. 2011) seems to support this (c.f., 0.96 in Wrayl6- 
423). 

In Figure 14a, we present the observed SED plots (blue dia¬ 


monds and grey lines) and the modelled SED (red line). The dashed 
line is the modelled incident SED of the CSPN, where the flux den¬ 
sity in the F574M band is consistent with the observed value. In 
Figure 14b, we display a close-up of the SED in the wavelength 
range covered by the Spitzer/IRS spectrum. The large mass frac¬ 
tion of the graphite grains is due to the rising mid-IR flux den¬ 
sity toward ~25 pm. Because we did not consider the possibility of 
very large PAH molecules discussed in Section 3.6.3, the 15-20 pm 
PAH plateau is not reproduced by the model. The complex of the 
10-14/jm PAHs and broad 11 pm feature in Wrayl6-423 could be 
reproduced with PAHs although the fitting is not perfect. The ob¬ 
served band profile of the 10-14 pm PAHs and broad 11 ^m feature 
is nearly flat top around 12 pm whereas the predicted band pro¬ 
file is very similar to the observed one except for the small lack 
of emission indicated by the arrow. The lack of emission could be 
filled if we add amorphous carbon grains in the model, because the 
Q a bs(A) of amorphous carbon calculated from the optical data of 
Rouleau & Martin (1991) appears in a small bump around 12/mi. 


4 DISCUSSION 

4.1 Progenitor mass estimate by comparison to AGB 
nucleosynthesis models 

Because the HDS spectrum shows a signature of [WC] CSPN in the 
C iv AA 5801/12 A band, the CSPN could be He-rich. Assuming that 
the CSPN is in the midst of He-burning, comparing the estimated L, 
and Teg on He-burning post-AGB evolution tracks with the initial 
Z=0.008 of Vassiliadis & Wood (1994) indicates a progenitor mass 
of 1.5-2.0 Mo. 

The presence of the Civ/LI 5801/12 A band could be ev¬ 
idence that Wrayl6-423 experienced a very late thermal pulse 
(VLTP). Vassiliadis & Wood (1994) demonstrated the He-burning 
post-AGB track for initially 1.5 M 0 stars, with Z=0.008; these stars 
experienced a VLTP at log lo r eff ~4.9 (-79430 K). The subsequent 
evolution was dominated by He-shell burning. The difference from 
normal He-burning evolution tracks is that there is no upturn lu¬ 
minosity, and the evolution after VLTP is very similar to the H- 
burning tracks. If Wray 16-423 experienced a VLTP, then the pro¬ 
genitor mass is estimated to be >1.5 M 0 . 

Next, we verify whether AGB nucleosynthesis models can ex¬ 
plain observed elemental abundances in an initial -1.5-2.0 Mo star. 
In Table 8, we compare the AGB model results of Karakas (2010) 
with our observation results. The 13 C pocket mass is not included. 
Besides He/C/N/O/Ne/S/Fe abundances, we list the ejected mass at 
the last thermal pulse and the core-mass of the CSPN in the last 
two lines of Table 8. Here, we assume that the estimated gas mass 
through our SED model is nearly equal to the ejected mass at the 
last thermal pulse. The gas mass estimated in our SED models de¬ 
pends on the outer nebula radius and the de-reddened H/3 flux. So 
far, the gas mass estimates in Wray 16-423 range from 0.248 M 0 
(Dudziak et al. 2000) to O.49±O.1OM 0 (see section 3.7.2). Among 
these model predictions, the 1.9M 0 model provides an excellent 
fit to the observed elemental abundances and the ejected mass. The 
consistency between the observed and the predicted Ne abundances 
indicates that the l3 C pocket may not form in Wrayl6-423. Even 
if we adopt the expected CEL C abundance of 8.70±0.14, as de¬ 
scribed in section 3.4, our conclusion does not change consider¬ 
ably. 
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Figure 14. (upper panel) Predicted SED of the PN Wrayl6-423 by Cloudy modeling (the red solid line and the dashed line). The red solid line is the emergent 
model spectrum of the CSPN plus the nebula, and the dashed line is the spectrum of the CSPN. The grey lines are the observed spectra taken by the Subaru/HDS 
and the Spitzer/ IRS. The blue diamonds are the photometry data of the CSPN plus the nebula from the GALEX (Fuv and Nuv bands), HST (F547M), 2MASS 
(near-IR JHKs), WISE (11.6 and 22.1 pm), and .4 KARI/FIS 65 pm. The green star is the photometry of the CSPN measured by the //.ST/(I : 547 VI ) image. 
(lower panel ) Closed-up plot for mid-IR wavelengths. The spectral resolving power ( R) of the Cloudy model spectrum was constant (100) to match that of the 
Spitzer/IRS. The lack of emission, indicated by the arrow, could be filled by adding amorphous carbon grains in the model. See the text for details. 




4.2 Expected solid-phase Mg, S, and Fe 

The difference between the observed S and Fe abundances and the 
AGB model predicted abundances could constrain how much mass 
of each element could be locked by the dust grains, based on the 
assumption that elemental abundances in the nebula are consistent 
with the predicted values in the case of the 1.9M o /Z=0.004 AGB 
model. The density fraction of an element X in the solid phase, 
n S oiid(X)/n(X), can be roughly estimated using the following: 


^solid(X) _ J _ 10 i ogl0 (g2) obs _iog 10 (g2) modd (g) 

n(X) 

where log 10 (/t(X)//!(H)) obs and log l0 (n(X)/n(H)) mO dei are the ob¬ 
served and the AGB star model predicted elemental abundances of 
the element X. Using the Equation (8), ~92-95 % of the Fe atoms 
and ~46-53 % of the S atoms in the nebula could be in the solid 
phase. 
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Table 7. Comparison between SED modeling and observations. 


Ion 

dlab 

/(CLOUDY) 

(/(H6)=100) 

/(Obs) 
(7(H/3)= 100.0) 

Ion 

-flab 

/(CLOUDY) 
(7(H/3)= 100.0) 

/(Obs) 
(7(H/3)= 100.0) 

[On] 

3726 A 

19.79 

19.52 

He i 

5876 A 

16.22 

15.20 

[On] 

3729 A 

14.17 

10.72 

[Sm] 

6312 A 

1.19 

1.42 

[Ne in] 

3869 A 

87.36 

85.71 

[Nn] 

6548 A 

5.45 

5.50 

[Ne in] 

3968 A 

26.33 

26.15 

Ha 

6563 A 

284.8 

285.0 

[S ii] 

4069 A 

0.63 

1.11 

[Nn] 

65 84 A 

16.08 

16.20 

Cm 

4187 A 

0.17 

0.15 

He i 

6678 A 

4.19 

3.98 

Cn 

4267 A 

0.45 

0.47 

[Sn] 

6716 A 

1.82 

1.58 

Hi5 

4341 A 

46.69 

46.90 

[Sn] 

6731 A 

2.26 

2.64 

[On.] 

4363 A 

14.59 

13.39 

[Arm] 

7135 A 

7.77 

10.83 

He i 

4388 A 

0.63 

0.57 

He i 

7281 A 

1.06 

1.01 

He i 

4471 A 

5.20 

5.10 

[On] 

7323 A 

1.56 

2.79 

[Fe hi] 

4701 A 

0.013 

0.022 

[On] 

7332 A 

1.24 

2.15 

[Ariv] 

4711 A 

2.81 

1.86 

[Arm] 

7751 A 

1.88 

2.71 

[Ar iv] 

4740 A 

3.72 

2.71 

[Cl iv] 

8047 A 

0.34 

0.56 

Cm 

4649 A 

0.48 

0.40 

[Sm] 

9069 A 

4.65 

7.99 

Civ 

4659 A 

0.045 

0.103 

[Siv] 

10.51 fj m 

55.07 

27.80 

Hen 

4686 A 

10.02 

11.14 

[Ne n] 

12.80/im 

0.21 

1.07 

H/3 

4863 A 

100.0 

100.0 

[Ne in] 

15.55/nn 

49.05 

59.76 

[Fern] 

4881 A 

0.032 

0.032 

[Sm] 

18.67 /im 

7.65 

9.92 

He i 

4922 A 

1.37 

1.27 

[Oiv] 

25.88 /im 

16.00 

11.92 

[O in] 

4931 A 

0.14 

0.14 

[Sm] 

33.47 /im 

3.99 

4.31 

[On.] 

4959 A 

336.3 

362.4 

[Nem] 

36.01 /im 

4.03 

3.87 

[On.] 

5007 A 

1012.1 

1062.1 

IRS-1 

6.20/im 

12.36 

13.466 

[Arm] 

5192 A 

0.12 

0.13 

IRS-2 

6.90//m 

11.29 

13.96 

[Fe in] 

5271 A 

0.064 

0.040 

IRS-3 

7.90/nn 

18.95 

18.22 

[Cl iv] 

5324 A 

0.012 

0.021 

IRS-4 

11.30/im 

19.45 

19.28 

[Cl hi] 

5518 A 

0.28 

0.25 

IRS-5 

12.30/im 

19.47 

21.30 

[Cl m] 

5538 A 

0.31 

0.37 

IRS-6 

13.50/im 

40.31 

40.70 

[Nn] 

5755 A 

0.45 

0.40 

IRS-7 

14.50/im 

21.20 

20.52 

Band 

A c 

7 v (CLOUDY) 

7 v (Obs) 

X 


l°g io( x / h )cloudy 

log,o(X/H)o bs 



(mJy) 

(mJy) 



+ 12 

+ 12 

GALEX -Fuv 

1516 A 

5.61 

2.58 

He 


11.02 

11.02 

GALEX -Nuv 

2267 A 

3.96 

3.07 

C 


8.86 

8.92 

WFPC2/F547M 

5484 A 

1.35 

1.44 

N 


7.69 

7.69 

2MASS-7 

1.235 /im 

2.49 

2.13 

O 


8.22 

8.31 

2MASS-/7 

1.662/im 

1.60 

1.15 

Ne 


7.50 

7.61 

2MASS-/G 

2.159/im 

2.47 

1.68 

S 


6.40 

6.40 

WISE- B1 

3.353/im 

1.90 

1.32 

Cl 


4.71 

4.74 

WISE-B2 

4.603 fi m 

2.87 

1.69 

Ar 


5.92 

6.00 

WISE- B3 

11.56/im 

23.79 

23.68 

Fe 


5.57 

5.56 

AKARI/F1S 

65 /im 

88.30 

<88.33 

^(reduced x 2 ) 



48.14(0.63) 


Note - The values for the IRS-1, 2, 3, 4, 5, 6, and 7 bands correspond to integrated fluxes between the following wavelength ranges, 5.9-6.9gm, 6.4-7.4//m. 
7.4-8.4/im, 10.8-11.8/im, 11.8-12.8/im, 12.5-14.5/nn, and 14-15/an, respectively. 


Table 8. Comparison to AGB nucleosynthesis models of Karakas (2010). 
The C abundance in Wrayl6-423 is the RL C value. See text regarding the 
estimate of the Mg abundance in Section 4.2. The core-mass of Wrayl6-423 
was estimated from Vassiliadis & Wood (1994). 




AGB models 

Present 

work 

1.5 M 0 
Z=0.004 

1.75 Mq 
Z=0.004 

1.90 M 0 
Z=0.004 

He 

10.97 

10.98 

10.99 

11.02+0.04 

C 

8.46 

8.80 

8.94 

8.92+0.08 

N 

7.65 

7.69 

7.71 

7.69+0.04 

O 

8.23 

8.25 

8.25 

8.31+0.01 

Ne 

7.42 

7.51 

7.61 

7.61+0.03 

Mg 

6.88 

6.88 

6.89 

6.19+0.30 

S 

6.69 

6.70 

6.70 

6.40+0.03 

Fe 

6.78 

6.78 

6.78 

5.56+0.07 

Ejected mass 

0.04 M 0 

0.15 M 0 

0.44 M 0 

0.59+0.10 M e 

Core-mass 

0.62 M 0 

0.63 M G 

0.64 M 0 

0.62-0.67 M q 


The ionised Mg abundances in PNe have been calculated 
using the forbidden lines at UV wavelengths and recombination 
Mg n ,14481 A. We were unable to detect optical Mgn lines. In¬ 
stead, we attempt to estimate the gas-phase Mg abundance using 
Mg° and N° abundances based on the assumption that all N atoms 
exist as a gas. This assumption is possible, because the observed N 
abundance (7.69±0.04) is consistent with AGB model predictions 
(7.71). Under the assumption that Mg-containing grains did not 
form, the Mg/N ratio would be 0.151 from the 1.9 M 0 AGB model; 
thus, nearly identical to Mg°/N°. In this case (i.e., Mg-atoms are 
not depleted by dust grains), the Mg°/H + should be 0.151xN°/H + = 
8.98(-8)±4.92(-9). Because the actual Mg°/H+ is 1.81(-8)±5.27(- 
9), about 14-26 % of the Mg-atoms likely exist as a gas and the re¬ 
maining 74-86 % are in the solid phase. Using the modeled Mg/H 
of 6.89 dex (7.76(-6) in linear scale), the gas-phase Mg/H would 
be 6.19±0.30 dex (1.53(-6)±4.63(-7) in linear scale). Both the gas 
and solid phase Mg estimates are optimistic and depend on the 
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Table 9. Required MgS and Fe50 grain masses for the 30 p m feature. The 
maximum Fe50S (Mgo. 5 Feo. 5 S) and MgS mass could form 7.65(—5) M 0 and 
5.95(-5) Mo, respectively. 


Models 

Trf(MgS) 

7T(Fe50S) 

m d ( MgS) 

»V(Fe50S) 


(K) 

(K) 

(Mo) 

(Mo) 

Case 1 

133.8±0.9 

117.9± 1.2 

5.5(-7)±1.0(-7) 

1. 8 (—6)±3.2(—7) 

Case 2 

// 

n 

6 . 0 (— 6 )± 1.1 (— 6 ) 

1.9(—5)±3.6(— 6 ) 


AGB model prediction. UV Mg forbidden lines in HST/STIS spec¬ 
tra would provide a more exact value. 

Using the estimated fractions of the solid phase Mg, S, and 
Fe atoms, we estimate the solid-phase masses m so ijd(X) using the 
following equation: 

m. n(X) n, n iid(X) 

iWX)=/r(X)- * -^ '^7^’ ( 9 ) 

1 + 4 n(He)/n(H) obs «(H) model n(X) 

where p(X) is the atomic mass, and u(He)/«(H) obs is the He 
number density (1.06(—1), See Table 4). 

In summary, ~92-95% of the Fe atoms, ~74-86% of the Mg 
and ~46-53% of the S atoms in the nebula could be in the solid 
phase, corresponding to masses of respectively 1.3±0.2, 0.7±0.2 
and 0.33±0.01 xl(T 4 M 0 . 


4.3 Can MgS or Fe50S grains explain the 30 gm feature? 

In largely Fe-depleted C-rich PNe, such as Wray 16-423, estimated 
carbon dust masses may make up a portion of the dust grains in the 
nebula, not including the broad 16-24pm and 30 pm features. As 
discussed in Section 3.6.4, Mg-Fe-S grains may form in the dusty 
nebula of Wray 16-423, and could be a carrier for the broad 30 pm 
feature. Mg-Fe-S grain masses are controlled by the solid-phase 
S-atorn mass. Adopting the solid-phase Mg, S, and Fe masses as 
we estimated above, the possible maximum Fe50S (Mgo. 5 Feo. 5 S) 
and MgS masses could form 7.65(—5) M 0 and 5.95(-5) M 0 , respec¬ 
tively. 

As such, can MgS or Fe50S grains explain the 30pm fea¬ 
ture? To answer this question, we estimate their masses using the 
Plank function with a single dust temperature method described in 
Hony et al. (2003), who explained that the dust mass m d in the op¬ 
tically thin nebula can be estimated using the following equation: 

4 r r\2 I r k Q»bs(A 0) /rr , lm 

m,i = - F MjM D p y --- B A (T d ) dAJ , (10) 

where F Ao j t is the integrated flux in the wavelength range from 
To to Ai . Here, we adopt the central wavelength and FWHM listed 
in Table Cl of Appendix C. T 0 and A t are 28pm and 36pm, re¬ 
spectively, D is the distance (24.8 kpc), and p is the grain den¬ 
sity, 3 g cm ' 3 for MgS and 5 g cnT 3 for Fe50S. When the grain 
radius is small in comparison with the wavelength T, the absorp¬ 
tion profile becomes independent of a (Pegourie 1988). We calcu¬ 
lated the Q a bs(A)/a values of the CDE MgS and Fe50S by follow¬ 
ing Equation (4) of Fabian et al. (2001). The optical data of MgS 
(Mgo. 9 Feo. 1 S, actually) and Fe50S (Mgo. 5 Feo. 5 S) were taken from 
Begemann et al. (1994). 

In the models, we use an average dust temperature for both 
MgS and Fe50S derived by Plank function fitting with a single 
temperature. For F A(h A ,, we consider two cases: (Case 1) the un¬ 
derlying continuum of the 30 pm feature is largely from graphite, 
as demonstrated in the Cloudy model, and (Case 2) instead of 


graphite, the contribution comes from amorphous carbon, whose 
absorption efficiency drops sharply for longer wavelengths, com¬ 
pared with graphite; thus, MgS or Fe50S may be the main con¬ 
tributors to this band radiation. For Case 1, we calculated F, 
of 4.88(-13) erg s' 1 cm' 2 using the residual spectrum generated 
by subtracting the Cloudy synthesised spectrum from the observed 
Spitzer /IRS spectrum. For Case 2, we calculated F Ao , Ai of 5.33(—12) 
erg s' 1 cm' 2 from the observed Spitzer/IRS spectrum (the Cloudy 
spectrum was not subtracted). 

The required masses to produce the band flux of the 30 pm 
feature are listed in the last column of Table 9. If the graphite grains 
exist in the nebula (Case 1), then the required MgS or Fe50S mass 
is less than the possible maximum MgS and Fe50S masses. Even if 
not graphite but amorphous carbon are the main dust grains in the 
nebula, MgS and Fe50S grains can reproduce the intensity peak of 
the observed broad 30pm feature (Case 2). 

At present, we hold onto the possibility that MgS and Fe50S 
are carriers of the 30pm feature in the case of Wrayl6-423. By 
adding the 1.8(-6)±3.2(-7) M Q of Fe50S, we could estimate the 
maximum dust mass to be 1.6(-2)M 0 . However, we understand 
that these estimates are optimistic, given that we did not consider 
S-atom-based molecules, such as CS. If any sulfide molecules form 
in the nebula, then it would be more difficult to explain the 30 pm 
feature with Fe50S and MgS. If it is possible to estimate the mass of 
the S-atom existing as the molecule with sub-mm spectra and the 
optical characteristics data of MgS, Fe50S and the new graphite 
data by Jiang et al. (2014) are available over a wider wavelength 
range, then it may be possible to resolve which grains carry the 
30pm feature. 


5 SUMMARY 

We performed a detailed chemical abundance and dust analysis of 
Sgr dSph PN Wrayl6-423. Numerous gas emission lines were de¬ 
tected, from which the abundances of 11 elements were calculated. 
The [Kr/H] and [Fe/H] abundances were comparable to Sgr dSph 
PN Hen2-436, which has a similar metallicity. The [Kr/O] abun¬ 
dance was consistent with the value expected by the [Kr/0]-C/0 
relation found amongst Galactic PNe. The extremely small [Fe/H] 
does not reflect the metallicity, but it implies that most of the Fe 
atoms would be in the solid phase, given the [Ar/H] abundance. The 
Spitzer !IRS spectrum displays broad 16-24pm and 30 pm features, 
as well as PAH bands at 6-9 pm (Class B) and 10-14 pm (Class a). 
An unidentified broad feature in ~ 16-24 pm is probably unrelated 
to FeS, amorphous silicate, and PAH. The progenitor was estimated 
to be a ~1.5-2.0 Mq star. The observed elemental abundances and 
the derived gas mass were in good agreement with the AGB nu¬ 
cleosynthesis model for an initial mass of 1.90 M 0 and a Z=0.004 
star. We found that about 80 %, 50 %, and 90 % of the Mg, S, and 
Fe atoms were likely in the solid phase, respectively. Using the de¬ 
rived solid-phase Mg, S, and Fe masses, we estimated the possible 
maximum MgS and Fe50S grain masses. We tested whether MgS 
and Fe50S grains can produce the intensity peak of the observed 
broad 30 pm band with these masses. Depending on the total sul¬ 
fide molecule masses, at present, MgS and Fe50S could be carriers 
of the 30pm feature in Wray 16-423. 
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Figure Al. Single 3680-5390 A HDS spectrum of Wray 16-423. The observed wavelength is corrected by the heliocentric radial velocity of +133.12 km s 1 . 
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Table Al. Detected lines and identifications in the HDS/FEROS spectra. BC in component 2 of C iv /15801/11 A means the broad component shown in Fig. 6. 
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2.718 

0.038 

6734.25 

[Su] 

6730.81 

2 

-0.320 

0.894 

0.026 

3870.72 

[Ne ml 

3869.06 

2 

0.291 

23.243 

1.367 

4771.41 

On 

4769.25 

1 

0.026 

0.017 

0.004 




T 


2.637 

0.071 




T 


85.713 

3.942 

4804.70 

Neu 

4802.58 

1 

0.016 

0.018 

0.003 

6747.42 

C in 

6744.39 

1 

-0.322 

0.027 

0.005 

3873.51 

He i 

3871.79 

1 

0.290 

0.101 

0.012 

4861.29 

Hen 

4859.32 

1 

0.001 

0.269 

0.095 

6783.13 

Cn 

6779.94 

1 

-0.326 

0.026 

0.004 

3890.10 

H8 

3889.05 

1 

0.287 

1.400 

0.123 

4861.74 

Hen 

4859.32 

2 

0.000 

0.430 

0.096 

6798.08 

IK iv] 

6795.10 

1 

-0.328 

0.030 

0.004 

3890.63 

H8 

3889.05 

2 

0.286 

17.390 

0.772 




T 


0.699 

0.135 

6830.07 

[Kr in] 

6826.85 

1 

-0.333 

0.030 

0.007 




T 


18.790 

0.834 

4863.50 

H4 

4861.32 

1 

0.000 

100.000 

0.882 

6894.08 

Hen 

6890.90 

1 

-0.341 

0.117 

0.011 

3922.37 

Cn 

3920.68 

1 

0.279 

0.034 

0.006 

4875.67 

N m 

4873.60 

1 

-0.003 

0.015 

0.006 

6992.04 

Cai 

6988.95 

1 

-0.354 

0.023 

0.003 

3925.24 

Hen 

3923.48 

1 

0.278 

0.097 

0.009 

4883.18 

[Fe in] 

4881.00 

1 

-0.005 

0.032 

0.002 

7009.13 

[Arv] 

7005.40 

1 

-0.356 

0.030 

0.003 

3928.30 

Hei 

3926.54 

1 

0.277 

0.116 

0.008 

4915.78 

Oi 

4913.60 

1 

-0.014 

0.040 

0.007 

7065.38 

He i 

7062.28 

1 

-0.364 

0.017 

0.003 

3966.48 

He i 

3964.73 

1 

0.267 

0.786 

0.034 

4924.11 

He i 

4921.93 

1 

-0.016 

1.270 

0.040 

7068.11 

He i 

7065.18 

1 

-0.364 

3.198 

0.148 

3969.10 

[Ne in] 

3967.79 

1 

0.267 

13.079 

0.884 

4926.86 

On 

4924.53 

1 

-0.017 

0.011 

0.004 

7068.59 

He i 

7065.18 

2 

-0.364 

4.287 

0.177 

3969.40 

[Ne in] 

3967.79 

2 

0.267 

13.068 

0.880 

4933.33 

[O in] 

4931.23 

1 

-0.019 

0.087 

0.003 




T 


7.484 

0.277 




T 


26.147 

1.438 

4933.76 

[O in] 

4931.23 

2 

-0.019 

0.031 

0.003 

7138.71 

[Ar in] 

7135.80 

1 

-0.374 

7.116 

0.159 

3971.83 

H7 

3970.07 

1 

0.266 

14.345 

0.608 




T 


0.118 

0.005 

7139.36 

[Ar in] 

7135.80 

2 

-0.374 

3.710 

0.093 

4011.03 

He i 

4009.26 

1 

0.256 

0.192 

0.010 

4955.35 

Nn 

4952.82 

1 

-0.024 

0.067 

0.034 




T 


10.826 

0.229 

4025.78 

He i 

4023.98 

1 

0.252 

0.033 

0.008 

4961.04 

[O in] 

4958.91 

1 

-0.026 

266.598 

3.656 

7180.68 

Hen 

7177.52 

1 

-0.379 

0.166 

0.028 

4027.87 

He i 

4026.18 

1 

0.251 

1.049 

0.055 

4961.46 

[O in] 

4958.91 

2 

-0.026 

95.810 

2.772 

7240.00 

Cai 

7236.81 

1 

-0.387 

0.494 

0.091 

4028.17 

He i 

4026.18 

2 

0.251 

1.035 

0.042 




T 


362.407 

4.612 

7284.59 

He i 

7281.35 

1 

-0.393 

1.011 

0.032 




T 


2.084 

0.087 

5008.93 

[O in] 

5006.84 

1 

-0.038 

581.111 

6.375 

7321.86 

[Ou] 

7318.92 

1 

-0.398 

0.546 

0.028 

4070.16 

[Sn] 

4068.60 

1 

0.239 

0.732 

0.028 

5009.34 

[O in] 

5006.84 

2 

-0.038 

481.000 

5.876 

7322.85 

[Ou] 

7318.92 

2 

-0.398 

1.292 

0.038 

4070.71 

[Sn] 

4068.60 

2 

0.239 

0.375 

0.017 




T 


1062.111 

8.960 

7323.67 

[Ou] 

7318.92 

3 

-0.398 

0.951 

0.037 




T 


1.107 

0.042 

5017.90 

He i 

5015.68 

1 

-0.040 

2.233 

0.028 




T 


2.789 

0.074 

4071.99 

On 

4069.89 

1 

0.239 

0.245 

0.014 

5049.96 

He i 

5047.74 

1 

-0.048 

0.239 

0.010 

7332.46 

[Ou] 

7329.66 

1 

-0.400 

0.754 

0.030 

4074.02 

On 

4072.15 

1 

0.238 

0.081 

0.009 

5133.34 

C hi 

5130.86 

1 

-0.067 

0.080 

0.006 

7333.44 

[Ou] 

7329.66 

2 

-0.400 

0.928 

0.039 

4077.89 

[Sn] 

4076.35 

1 

0.237 

0.284 

0.013 

5160.89 

[Fen] 

5158.78 

1 

-0.073 

0.010 

0.002 

7334.33 

[Ou] 

7329.66 

3 

-0.400 

0.467 

0.031 

4078.46 

[Sn] 

4076.35 

2 

0.237 

0.076 

0.005 

5182.38 

C hi 

5179.90 

1 

-0.078 

0.021 

0.005 




T 


2.150 

0.068 




T 


0.361 

0.015 

5194.02 

[Ar in] 

5191.82 

1 

-0.081 

0.126 

0.006 

7533.84 

Cn 

7530.57 

1 

-0.426 

0.224 

0.018 

4099.06 

N in 

4097.35 

1 

0.231 

0.355 

0.016 

5199.84 

[Ni] 

5197.90 

1 

-0.082 

0.054 

0.002 

7750.93 

Fei 

7747.60 

1 

-0.455 

0.140 

0.015 

4099.35 

N in 

4097.35 

2 

0.231 

0.157 

0.010 

5200.57 

[Ni] 

5197.90 

2 

-0.082 

0.027 

0.003 

7754.24 

[Ar in] 

7751.10 

1 

-0.455 

1.535 

0.051 




T 


0.513 

0.022 




T 


0.080 

0.003 

7754.94 

[Ar in] 

7751.10 

2 

-0.455 

1.172 

0.046 

4101.87 

Hen 

4100.04 

1 

0.230 

0.117 

0.043 

5202.12 

[Ni] 

5200.26 

1 

-0.083 

0.023 

0.002 




T 


2.707 

0.081 

4103.42 

H6 

4101.73 

1 

0.230 

12.501 

0.454 

5202.73 

[Ni] 

5200.26 

2 

-0.083 

0.033 

0.004 

7806.18 

[Viu] 

7802.76 

1 

-0.462 

0.077 

0.023 

4103.72 

H6 

4101.73 

2 

0.229 

10.173 

0.394 




T 


0.057 

0.004 

7819.78 

He i 

7816.14 

1 

-0.464 

0.163 

0.035 




T 


22.675 

0.802 

5243.20 

Nn 

5240.86 

1 

-0.091 

0.018 

0.004 

8011.56 

Fern 

8007.90 

1 

-0.488 

0.148 

0.031 

4105.19 

N in 

4103.39 

1 

0.229 

0.183 

0.016 

5272.75 

[Fe ui] 

5270.40 

1 

-0.098 

0.040 

0.005 

8025.93 

He i 

8021.30 

1 

-0.490 

0.141 

0.018 

4121.13 

On 

4119.22 

1 

0.224 

0.036 

0.007 

5325.45 

[Cl IV] 

5322.99 

1 

-0.108 

0.021 

0.004 

8049.26 

[Cl IV] 

8046.30 

1 

-0.492 

0.557 

0.026 

4122.66 

He i 

4120.81 

1 

0.224 

0.271 

0.011 

5344.85 

Cn 

5342.43 

1 

-0.112 

0.042 

0.003 

8240.55 

Hen 

8236.79 

1 

-0.515 

0.460 

0.038 

4130.60 

Nn 

4128.66 

1 

0.222 

0.034 

0.006 

5348.29 

[Kr iv] 

5346.02 

1 

-0.113 

0.093 

0.008 

8302.46 

P28 

8298.83 

1 

-0.522 

0.160 

0.020 

4145.59 

He i 

4143.76 

1 

0.217 

0.342 

0.013 

5377.63 

Neu 

5375.14 

1 

-0.119 

0.014 

0.002 

8310.10 

P27 

8306.11 

1 

-0.523 

0.204 

0.026 

4158.28 

Nn 

4156.36 

1 

0.213 

0.029 

0.006 

5413.96 

Hen 

5411.52 

1 

-0.126 

1.036 

0.023 

8317.78 

P26 

8314.26 

1 

-0.524 

0.235 

0.028 

4164.82 

C in 

4162.88 

1 

0.211 

0.020 

0.006 

5440.51 

Cui 

5438.08 

1 

-0.131 

0.022 

0.005 

8327.03 

P25 

8323.42 

1 

-0.525 

0.207 

0.020 

4170.84 

He i 

4168.97 

1 

0.210 

0.048 

0.005 

5520.13 

[Cl in] 

5517.72 

1 

-0.145 

0.251 

0.007 

8337.53 

P24 

8333.78 

1 

-0.526 

0.261 

0.025 

4182.98 

Ar ii 

4181.23 

1 

0.206 

0.016 

0.005 

5540.28 

[Cl in] 

5537.89 

1 

-0.149 

0.365 

0.009 

8349.30 

P23 

8345.55 

I 

-0.527 

0.205 

0.018 

4188.74 

C in 

4186.90 

1 

0.204 

0.150 

0.009 

5579.42 

[i in] 

5576.94 

1 

-0.155 

0.015 

0.004 

8362.67 

P22 

8359.00 

1 

-0.529 

0.297 

0.025 

4188.90 

N in 

4187.06 

1 

0.204 

0.314 

0.059 

5756.82 

[Nn] 

5754.64 

1 

-0.185 

0.293 

0.010 

8365.14 

He i 

8361.73 

1 

-0.529 

0.136 

0.020 

4197.62 

N in 

4195.74 

1 

0.202 

0.015 

0.007 

5757.56 

[Nn] 

5754.64 

2 

-0.185 

0.106 

0.007 

8371.16 

Piv 

8367.50 

1 

-0.530 

0.078 

0.018 

4201.73 

Hen 

4199.83 

1 

0.200 

0.236 

0.010 




T 


0.399 

0.013 

8378.14 

P21 

8374.48 

1 

-0.531 

0.294 

0.018 

4229.54 

[Fev] 

4227.19 

1 

0.192 

0.040 

0.005 

5804.05 

C IV 

5801.34 

1 

-0.192 

0.053 

0.004 

8396.10 

P20 

8392.40 

1 

-0.533 

0.297 

0.017 

4269.06 

Cn 

4267.18 

1 

0.180 

0.469 

0.015 

5808.55 

C iv 

BC 

2 

-0.193 

0.660 

0.056 

8417.00 

P19 

8413.32 

1 

-0.535 

0.302 

0.019 

4305.64 

On 

4303.82 

1 

0.168 

0.018 

0.006 

5814.59 

C iv 

5012.00 

3 

-0.193 

0.066 

0.004 

8441.71 

P18 

8437.95 

1 

-0.537 

0.368 

0.015 

4327.74 

On 

4325.76 

1 

0.161 

0.037 

0.008 




T 


0.779 

0.056 

8470.98 

P17 

8467.25 

1 

-0.541 

0.428 

0.018 

4340.64 

Hen 

4338.67 

1 

0.157 

0.341 

0.234 

5870.36 

[Kr iv] 

5867.74 

1 

-0.202 

0.111 

0.005 

8506.15 

P16 

8502.48 

1 

-0.544 

0.535 

0.033 

4342.39 

H5 

4340.46 

1 

0.157 

46.900 

0.724 

5878.10 

He i 

5875.60 

1 

-0.203 

9.704 

0.227 

8549.06 

P15 

8545.38 

1 

-0.549 

0.659 

0.048 

4351.39 

On 

4349.43 

1 

0.154 

0.022 

0.004 

5878.59 

He i 

5875.60 

2 

-0.203 

5.494 

0.172 

8602.21 

P14 

8598.39 

1 

-0.554 

0.738 

0.032 

4365.13 

[O in] 

4363.21 

1 

0.149 

13.397 

0.208 




T 


15.198 

0.331 

8668.84 

P13 

8665.02 

1 

-0.560 

0.937 

0.036 

4381.09 

N in 

4379.20 

1 

0.144 

0.039 

0.005 

5899.35 

Hen 

5896.78 

1 

-0.206 

0.025 

0.003 

8754.33 

P12 

8750.47 

1 

-0.568 

1.172 

0.040 

4389.75 

He i 

4387.93 

1 

0.142 

0.302 

0.011 

5955.65 

Hen 

5952.93 

1 

-0.214 

0.024 

0.004 

9018.92 

P10 

9014.91 

1 

-0.590 

1.874 

0.061 

4390.06 

He i 

4387.93 

2 

0.142 

0.269 

0.012 

6007.32 

Hen 

6004.72 

1 

-0.222 

0.022 

0.005 

9071.93 

[S in] 

9068.60 

1 

-0.594 

0.759 

0.071 




T 


0.570 

0.018 

6039.29 

Hen 

6036.78 

1 

-0.226 

0.018 

0.003 

9072.70 

[S in] 

9068.60 

2 

-0.594 

5.244 

0.215 

4394.14 

Neu 

4391.99 

1 

0.140 

0.019 

0.005 

6076.85 

Hen 

6074.19 

1 

-0.232 

0.017 

0.003 

9073.07 

[S in] 

9068.60 

3 

-0.594 

1.983 

0.613 

4439.48 

He i 

4437.55 

1 

0.126 

0.073 

0.005 

6104.44 

[K iv] 

6101.79 

1 

-0.235 

0.141 

0.004 




T 


7.985 

0.653 

4473.48 

He i 

4471.47 

1 

0.115 

5.167 

0.126 

6120.93 

Hen 

6118.26 

1 

-0.238 

0.038 

0.005 








4512.86 

[K iv] 

4510.92 

1 

0.103 

0.024 

0.003 

6173.48 

Hen 

6170.69 

1 

-0.245 

0.045 

0.007 
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Table A2. Line and band flux measurements in the Spitzer spectrum. 


4obs 

(pm) 

Ions/ 

PAH 

4lab 

(pm) 

H/l)” 

(erg s 1 cm 2 ) 

m 

umw 

(/(H8)=100) 

10.50 

LSivJ 

10.51 

5.42(—13)±3.63(—14) 

-0.959 

27.796±2.340 

11.28 

PAH 

11.30 

4.39(-14)±3.40(-15) 

-0.970 

2.245±0.208 

12.38 

Hi 

12.40 

4.44(-14)±5.68(-15) 

-0.980 

2.266±0.312 

12.80 

[Ne n] 

12.80 

2.10(—14)±2.45(—15) 

-0.983 

1.069±0.136 

15.56 

[Ne hi] 

15.56 

1. 17(-12)±2.25(-14) 

-0.985 

59.760±3.249 

18.73 

[Sm] 

18.71 

1,94(—13)±5.97(—15) 

-0.981 

9.917±0.589 

25.92 

[Oiv] 

25.89 

2.34(-13)±7.29(-15) 

-0.989 

11.921 ±0.711 

33.45 

[Sm] 

33.50 

8.49(-14)±2.31(-15) 

-0.993 

4.314±0.249 

36.02 

[Ne hi] 

36.00 

7.62(-14)±2.83(-15) 

-0.993 

3.874±0.244 


Note - F(H/3) is (1.53±0.02)xl(r 12 erg s" 1 cm~ 2 . 



Wavelength (/xm) Wavelength (/xm) 


Figure Cl. (left panels ) The 5.5-14/im continuum subtracted spectra, 
showing the 6-9 pm and 10-14/im PAH bands and the broad 11 pm fea¬ 
ture. Spitzer/IRS continuum subtracted spectrum of the Galactic PN K3-60, 
indicated by the grey line in each panel, is presented as a template for the 
16.4^/m PAH and 15-18/rm plateau. We convolved the SH and LH spec¬ 
tra with a Gaussian kernel for the SH and LH module spectra of K3-60 to 
smooth its spectral resolving power down to those of the SL and LL spectra, 
respectively. After that, we subtracted the local continuum and gas-emission 
lines from this spectrum, and then normalised to the intensity peak of the 
7.7 pm PAH band. (right panels) The 14-27 pm continuum subtracted spec¬ 
tra, showing the broad 16-24/^m feature. The spectrum of K3-60, indicated 
by the grey line, is superimposed on each spectrum of the sample PNe. The 
dashed lines show the wavelength range of the 15-18/im plateau + 16.4/im 
PAH. 


Table Bl. Ionic abundances from CELs. 


X m+ 

^lab 

/U lab )//(H/J) 

[7(H/3)=100] 

T c 

(K) 

n e 

(cm -3 ) 

X m+ /H + 

N° 

5197.90 A 

8.04(—2)±3.29(—3) 

10000 

1010 

6.11(-7)±2.50(-8) 


5200.26 A 

5.65(—2)±4.29(—3) 

10000 

1010 

5.71(—7)±4.33(—8) 

5.94(-7)±3.26(-8) 

N + 

5754.64 A 

3.99(—1)±1.25(—2) 

12090 

3540 

2.15(—6)±1.78(—7) 


6548.04 A 

5.50(0)±1.40(-1) 

12090 

3540 

2.08(-6)±1.03(-7) 


6583.46 A 

1.62(+1)±4.13(-1) 

12090 

3540 

2.07(-6)±1.02(-7) 

2.07(-6)±1.04(-7) 

0 + 

3726.03 A 

1.95(+1)±9.38(-1) 

12090 

3540 

8.33(—6)±6.83(—7) 


3728.81 A 

1.07(+1)±5.31(-1) 

12090 

3540 

8.39(—6)±7.08(—7) 


7320/7330 A 4.08(0)±2.57(-l) 

12090 

6870 

1.07(-5)±1.17(-6) 

8.63(-6)±7.49(-7) 

0 2+ 

4363.21 A 

1.34(+1)±2.08(-1) 

12380 

8820 

1.90(-4)±7.30(-6) 


4931.23 A 

1.18(—1)±4.52(—3) 

12380 

8820 

1.57(—4)±6.67(—6) 


4958.91 A 

3.62(+2)±4.61(0) 

12380 

8820 

1.88(—4)±4.22(—6) 


5006.84 A 

1.06(+3)±8.96(0) 

12380 

8820 

1.91(-4)±3.91(-6) 

1.90(-4)±4.02(-6) 

0 3+ 

25.89 pm 

1.19(+1)±7.11(-1) 

12800 

11320 

6.93(-6)±4.33(-7) 

Ne + 

12.81pm 

1.07(0)±1.36(-1) 

12070 

7570 

1.31(-6)±1.71(-7) 

Ne 2+ 

3869.06 A 

8.57(+l)±3.94(0) 

12800 

11320 

3.75(-5)±3.26(-6) 


3967.79 A 

2.61(+1)±1.44(0) 

12800 

11320 

3.80(—5)±3.49(—6) 


15.56pm 

5.98(+l)±3.25(0) 

12800 

11320 

3.76(-5)±2.10(-6) 


36.00 pm 

3.87(0)±2.44(-l) 

12800 

11320 

3.26(—5)±2.11(—6) 

3.7S(-5)±2.87(-6) 

Ne 3+ 

4724.17 A 

1.28(-2)±2.36(-3) 

12800 

11320 

2.23(—6)±5.34(—7) 


4725.64 A 

1.09(-2)±2.88(-3) 

12800 

11320 

1.98(—6)±6.07(—7) 

2.12(-6)±5.68(-7) 

Mg° 

4571.10A 

2.93(-l)±8.53(-2) 

10000 

1010 

1.81(-8)±S.27(-9) 

S + 

4068.60 A 

1.11 (0)±4.17(—2) 

12090 

3970 

1.25(—7)±8.07(—9) 


6716.44 A 

1.58(0)±5.01(-2) 

12090 

3970 

1,07(—7)± 1.13(—8) 


6730.81 A 

2.64(0)±7.14(-2) 

12090 

3970 

1.07(-7)±1.07(-8) 

l.ll(-7)±1.03(-8) 

S 2+ 

6312.10 A 

1.42(0)±3.44(-2) 

12070 

6860 

1.61(-6)±2.72(-7) 


9068.60 A 

7.99(0)±6.53(-l) 

12070 

6860 

1.43(-6)±1.68(-7) 


18.71pm 

9.92(0)±5.89(-l) 

12070 

6860 

1.44(-6)±1.04(-7) 


33.50 pm 

4.31 (0)±2.49(—1) 

12070 

6860 

1.44(-6)± 1.04(—7) 

1.44(-6)±1.36(-7) 

S 3+ 

10.51 pm 

2.78(+l)±2.34(0) 

12380 

8820 

9.57(-7)±8.06(-8) 

Cl 2+ 

5517.72 A 

2.51(—1)±6.81(—3) 

12070 

8280 

2.65(—8)±3.35(—9) 


5537.89 A 

3.65(—1)±8.60(—3) 

12070 

8280 

2.64(—8)±3.26(—9) 

2.65(-8)±3.30(-9) 

Cl 3+ 

5322.99 A 

2.08(-2)±4.28(-3) 

11800 

8820 

2.68(-8)±9.79(-9) 


8046.30 A 

5.57(—1)±2.65(—2) 

11800 

8820 

2.78(—8)±4.41(—9) 

2.78(-8)±4.60(-9) 

Ar 2+ 

5191.82A 

1.26(—1)±5.62(—3) 

12000 

7570 

6.70(—7)±6.46(—8) 


7135.80A 

1.08(+1)±2.29(-1) 

12000 

7570 

6.66(—7)±3.09(—8) 


7751.10A 

2.71 (0)±8.10(—2) 

12000 

7570 

6.94(—7)±3.55(—8) 

6.71(-7)±3.21(-8) 

Ar 34 

4711.37 A 

1.86(0)±2.64(-2) 

12800 

11320 

3.25(-7)±2.18(-8) 


4740.16 A 

2.72(0)±3.80(-2) 

12800 

11320 

3.28(—7)±2.10(—8) 

3.27(-7)±2.13(-8) 

Ar 4+ 

7005.40 A 

2.99(-2)±3.40(-3) 

12800 

11320 

3.11(-9)±3.87(-10) 

k 3+ 

6101.79 A 

1.41(-l)±4.42(-3) 

12800 

11320 

6.62(—9)±4.84(—10) 


6795.10 A 

3.02(-2)±3.78(-3) 

12800 

11320 

6.52(—9)±9.23(—10) 

6.60(-9)±S.61(-10) 

Fe 2+ 

4701.53 A 

2.21(—2)±2.72(—3) 

12090 

3540 

1.49(-8)±1.98(-9) 


4881.00 A 

3.17(—2)±2.07(—3) 

12090 

3540 

1,05(—8)±8.89(—10) 


5270.40 A 

4.00(-2)±5.45(-3) 

12090 

3540 

1.91 (—8)±2.74(—9) 

1.52(-8)±1.94(-9) 

Kr 24 

6826.85 A 

2.99(—2)±7.37(—3) 

12070 

7570 

1.81(—9)±4.81(—10) 

Ki J+ 

5867.74 A 

1.11(—1)±5.07(—3) 

12380 

8820 

1.71(—9)±8.18(—11) 


Note - Corrected recombination contributions for the [O n] AA 7320/30 A 
and for the [O m] A 4363 A lines. 
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Table B2. Ionic abundances from RLs. 



4lab 

Multi. 

A/llab) 

[/(H£)=100] 

X m+ /H + 

He + 

5875.60 A 

Vll 

1.52(+1) ± 3.31(—1) 

9.566—2) 

± 6.75(-3) 


4471.47 A 

V14 

5.17(0) ± 1.26(—1) 

9.68(-2) 

± 2 . 12 (- 2 ) 


6678.15 A 

V46 

3.97(0) ± 1.12(— 1) 

9.89(-2) 

± 6.35(-3) 


4921.93 A 

V48 

1.27(0) ±3.98(-2) 

9.55(-2) 

± 5.47(-3) 


4387.93 A 

V51 

5.706—1) ± 1.81 (—2) 

9.60(-2) 

9.63(-2) 

± 3.05(-3) 

± 9.40(-3) 

He 2+ 

4685.71 A 

V3.4 

l.U(+l) ± 1.57(-1) 

9.35(-3) 

± 2.16(-3) 

C 2+ 

6578.05 A 

V2 

2.11(-1) ± 7.76(-3) 

4.70(-4) 

± 1.63 (-4) 


4267.18 A 

V 6 

4.69(—1) ± 1.46(-2) 

4.73(-4) 

4.72(-4) 

± 1.29(-4) 

± 1.40(-4) 

C 3+ 

4647.42 A 

VI 

2.31(—1) ± 4.84(-3) 

3.65 (-4) 

± 7.48(-5) 


4650.25 A 

VI 

1.73(—1) ± 5.406—3) 

3.43(-4) 

± 7.07(-5) 


4186.90 A 

V18 

1.506-1) ± 8.676-3) 

2.69(-4) 

3.32(-4) 

± 6.86(-5) 

± 7.18(-5) 

C 4 + 

4658.20 A 

V 8 

1.036—1) ± 1.136-2) 

2.45(-5) 

± 6.75(-6) 

0 2+ 

4638.86 A 

VI 

4.15(—2) ± 4.79(-3) 

3.26(-4) 

± 8.15(—5) 


4649.13 A 

VI 

1.116-D ± 4.136-3) 

2.44(-4) 

± 5.43(-5) 


4661.63 A 

VI 

3.44(-2) ± 1.696—3) 

2.38(-4) 

± 5.43(-5) 


4676.23 A 

VI 

3.73(—2) ± 3.566-3) 

3.90(-4) 

± 9.56(-5) 


4349.43 A 

V2 

2.15(—2) ± 4.01(-3) 

2.32(-4) 

± 6.62(-5) 


4072.15 A 

V10 

8.08(-2) ± 8.75(—3) 

3.44(-4) 

± 8.50(-5) 


4119.22A 

V20 

3.6K-2) ± 7.37(-3) 

4.29(-4) 

± 1.27(—4) 


4924.53 A 

V28 

1.156-2) ± 4.156-3) 

2.66(-4) 

± 1.13(-4) 


4303.82 A 

V53a 

1.766-2) ± 5.83(—3) 

4.18(—4) 

3.12(-4) 

± 1.72(-4) 

± 8.19(-5) 


Table B3. Adopted ionisation correction factors (ICFs). 


X 

Line 

ICF(X) 

X/H 

He 

RL 

1 

He + +He 2+ 

C 

RL 

1 

c 2 + +c 3 + +c 4 

N 

CEL 

(^)cEL 

ICF(N)N + 

O 

CEL 

1 

o + +o 2 + +o 3+ 


RL 

(o^)cEL 

ICF(0)0 2+ 

Ne 

CEL 

1 

Ne + + Ne 2+ + Ne 3+ 

S 

CEL 

1 

S + + S 2+ + S 3+ 

Cl 

CEL 

1 

C1 2 + +C1 3+ 

Ar 

CEL 

1 

Ar 2+ + Ar 3+ 

K 

CEL 

(£0 

ICF(K)K 3+ 

Fe 

CEL 

(^)cEL 

ICF(Fe)Fe 2+ 

Kr 

CEL 

1 

Kr 2 + + Kr 3+ 
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Table Cl. Classification of the 6.2/7.7/8.6/im PAH and the 10-14/zm PAH bands and the measurements of the broad 11/16-24/30 /urn features. 


PNe 

6 .2//m 

PAH 

Class 

7.7 /dm 

PAH 

Class 

8.6 jum 

PAH 

Class 

10-14/im 

PAH 

Class 

11 p m 

Ac 

(pm) 

feature 

FWHM 

(pm) 

Ac 

(pm) 

m feature 

FWHM 

(pm) 

30 pm 

Ac 

(pm) 

feature 

FWHM 

(pm) 

LMC19 

8 


8 

a 

11.98±0.06 

1.85±0.17 

19.32±0.06 

6.09±0.25 

31.83±0.05 

5.58±0.22 

LMC25 

8 


8 

s 

11.46±0.01 

1.91 ±0.03 

18.47±0.03 

4.34±0.14 

29.82±0.02 

6.74±0.09 

LMC48 

81 


8 

6 

11.46±0.01 

1.82±0.02 

18.53±0.06 

4.15±0.11 

30.29±0.02 

7.08±0.06 

LMC78 

8 


8 

a 

11.87±0.02 

1.62±0.04 

18.81 ±0.05 

4.39±0.11 

31.43±0.05 

5.43±0.12 

Wray 16-423 

8 

8 

? 

a 

11,90±0.06 

1,76±0.11 

19.50±0.03 

6.47±0.08 

31.87±0.06 

5.57±0.13 


Table C2. The ratios of the 6-9 /dm PAH band, the broad 16-24 /dm feature, 
the broad 30 /im feature, and the 11 /dm SiC fluxes in the local continuum 
subtracted spectra to the integrated flux (F(IR)) between 5.75-36.8//m. 


PNe 

P(6-9pm PAH)/ F( 16-24pm)/ F(30pm)/ F(llpmSiC)/ 

F(IR) F(IR) F(IR) F(IR) 

LMC8 

1.96(^2) 

5.43C-2) 

2.45(-2) 

9.0K-2) 

LMC19 

7.28(-2) 

4.75(-2) 

3.13(—2) 


LMC25 

4.14(—2) 

2 .6K-2) 

5.52(—2) 

7.05(-2) 

LMC48 

5.19(-2) 

2.64(-2) 

4.85(-2) 

6 .88(-2) 

LMC78 

8.40(-2) 

2.55(-2) 

2.92(-2) 


LMC85 

5.58(-2) 

3.39(-2) 

2.41(—2) 

1-13(—1) 

LMC99 

1.07(—1) 

4.17(-2) 

1.29(—2) 


Wray 16-423 

1.77(—2) 

8.23(-2) 

3.33(-2) 
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